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of  the  two  disks.  Using  a  dyadic  Green’s  function  formulation  for  horizontally  stratified 
media  and  the  vector  Hankel  transform,  the  mixed  boundary  value  problem  is  reduced  to 
a  set  of  coupled  vector  integral  equations.  Employing  Galerkin’s  method  in  the  spectral 
domain,  the  complex  resonant  frequencies  are  calculated  and  convergence  of  the  results  is 
demonstrated.  It  is  shown  that  for  each  mode,  the  stacked  circular  microstrip  structure  has 
dual  resonant  frequencies  which  are  associated  with  the  two  coupled  constitutive  resonators 
of  the  structure  and  which  are  a  function  of  the  mutual  coupling  between  them.  This  mutal 
coupling  depends  on  the  geometrical  configuration  of  the  stacked  structure,  the  layered 
parameters,  and  the  disk  radii.  The  maximum  coupling  effect  occurs  where  the  real  parts 
of  the  resonant  frequencies  of  the  constitutive  resonators  are  approximately  equal,  where 
the  behavior  of  the  resonances  in  this  region  is  a  function  of  the  coupling.  The  dua1 
frequency  behavior  of  the  stacked  microstrip  structure,  easily  controlled  by  varying  the 
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parameters  of  layer  2  and  disk  radii  ratio,  given  fixed  parameters  for  layer  1  and  layer  3, 
may  be  used  to  broaden  the  bandwidth  or  provide  for  dual  frequency  use  of  the  antenna. 

A  rigorous  dyadic  Green’s  function  formulation  in  the  spectral  domain  is  used  to 
study  the  dispersion  characteristics  of  signal  strip  lines  in  the  presence  of  metallic  crossing 
strips.  A  set  of  coupled  vector  integral  equations  for  the  current  distribution  on  the 
conductors  is  derived.  Galerkin’s  method  is  then  applied  to  derive  the  matrix  eigenvalue 
equation  for  the  propagation  constant.  The  dispersion  properties  of  the  signal  lines  are 
studied  for  both  cases  of  finite  and  infinite  length  crossing  strips.  The  effects  of  the 
structure  dimensions  on  the  passband  and  stopband  characteristics  are  investigated.  For 
crossing  strips  of  finite  length,  the  stopband  is  mainly  affected  by  the  period,  the  crossing 
strip  length,  and  the  separation  between  the  signal  and  the  crossing  strips.  For  crossing 
strips  of  infinite  length  carrying  travelling  waves,  attenuation  along  the  signal  line  exists 
over  the  whole  frequency  range  of  operation. 

A  new  method  for  analyzing  frequency-dependent  transmission  line  systems  with 
nonlinear  terminations  is  presented.  The  generalized  scattering  matrix  formulation  is  used 
as  the  foundation  for  the  time  domain  iteration  scheme.  Compared  to  the  admittance 
matrix  approach  proposed  in  a  previous  paper,  it  has  the  advantage  of  shorter  impulse 
response  which  leads  to  smaller  computer  memory  requirement  and  faster  computation 
time.  Examples  of  a  microstrip  line  loaded  with  nonlinear  elements  are  given  to  illustrate 
the  efficiency  of  this  method. 
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RESONANT  FREQUENCIES 
OF  STACKED  CIRCULAR  MICROSTRIP  ANTENNAS 
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Abstract-  The  complex  resonant  frequencies  of  a  microstrip  antenna  consisting  of  two 
circular  microstrip  disks  in  a  stacked  configuration  are  investigated.  Using  a  dyadic  Green’s 
function  formulation,  a  rigorous  analysis  of  the  two  stacked  circular  disks  in  a  layered 
medium  is  performed.  The  mixed  boundary  value  problem  is  reduced  to  a  set  of  coupled 
vector  integral  equations  using  the  vector  Hankel  transform.  Employing  Galerkin’s  method 
in  the  spectral  domain,  the  complex  resonant  frequencies  are  calculated  as  a  function  of 
the  layered  parameters  and  the  ratio  of  the  radii  of  the  two  disks.  It  is  shown  that  for  each 
mode  the  stacked  circular  microstrip  structure  has  dual  resonant  frequencies  which  may 
provide  for  dual  frequency  or  wide  bandwidth  use. 
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1.  INTRODUCTION 
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Conventional  microstrip  antennas,  consisting  of  a  single  perfectly  conducting  patch  on 
a  grounded  dielectric  slab,  have  received  much  attention  in  recent  years  [l]-[2]  due  to  their 
many  advantanges,  including  low  profile  and  light  weight.  However,  due  to  their  resonant 
behavior,  their  use  is  severely  limited  in  that  they  radiate  efficiently  only  over  a  narrow 
band  of  frequencies,  with  bandwidths  typically  only  a  few  percent  [1].  Techniques  for 
increasing  the  bandwidth  have  included  stacking  a  number  of  microstrip  patches  in  mul¬ 
tilayer  configurations,  introducing  additional  resonances  in  the  frequency  range  of  interest 
and  achieving  bandwidths  up  to  approximately  10-20  percent. 

The  first  multilayered  microstrip  element  was  described  by  Oltman  [3], [4]  as  an  elec- 
tromagnetically  coupled  microstrip  dipole  where  a  printed  dipole  was  excited  by  an  open 
ended  microstrip  transmission  line  in  the  same  plane  as  the  dipole  or  in  the  layer  below  the 
dipole.  Hall  et  al.  [5]  stacked  rectangular  microstrip  patches  in  two-  and  three-layer  config¬ 
urations,  achieving  bandwidths  in  excess  of  16  times  that  of  alumina  substrate  microstrip 
antennas,  and  noted  that  the  stacked  configurations  allowed  for  simple  antenna/circuit 
integration.  Experimental  work  by  Sabban  [6]  and  Chen  et  al.  [7]  with  two-layer  stacked 
circular  microstrip  patches  produced  wider  bandwidths  and  higher  efficiencies  than  those 
obtained  with  conventional  single  patch  configurations.  Stacking  microstrip  patches  for 
dual  frequency  use  was  investigated  experimentally  for  circular  disks  by  Long  et  al.  [8] 
and  for  annular  rings  by  Dahele  et  al.  [9]. 

Rigorous  analyses  employing  Galerkin’s  method  have  been  applied  to  conventional 
single-layer  microstrip  antennas  of  circular,  rectangular,  annular  ring,  and  elliptic  geome¬ 
tries  [10]-[14j.  The  coupling  between  two  microstrip  antennas  on  a  single  substrate  has 
also  been  rigorously  treated  [15]-[17] .  However,  theoretical  analyses  of  stacked  microstrip 
patches  is  limited.  Printed  microstrip  dipoles  in  stacked  configurations  coupled  to  an  open 
ended  microstrip  transmission  line  have  been  studied  by  Katehi  et  al.  [18].  Tulintseff  et  al. 
[19]  used  the  method  of  moments  with  triangular  basis  functions  to  analyze  the  open  struc¬ 
ture  of  a  two-layer  circular  microstrip  antenna  excited  by  an  incident  plane  wave.  Kastner 
et  al.  [20]  described  a  spectral  domain  iterative  anlaysis  of  single-  and  double-layered 
microstrip  antennas  algorithm  to  compute  radiation  patterns. 

Considered  here  is  a  microstrip  antenna  consisting  of  two  circular  microstrip  disks  in 
a  stacked  configuration.  The  two  stacked  configurations  shown  in  Figures  1(a)  and  1(b), 
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denoted  configurations  A  and  5,  respectively,  are  investigated.  A  rigorous  analysis  of  the 
two  stacked  circular  disks  in  a  layered  medium  is  performed  using  a  dyadic  Green’s  function 
formulation.  Using  the  Vector  Hankel  transform,  the  mixed  boundary  value  problem  is 
reduced  to  a  set  of  coupled  vector  integral  equations  and  solved  by  employing  Galerkin’s 
method  in  the  spectral  domain.  The  current  distribution  on  each  disk  is  expanded  in 
terms  of  a  complete  set  of  transverse  magnetic  (TM)  and  transverse  electric  (TE)  modes 
of  a  cylindrical  resonant  cavity  with  magnetic  side  walls.  Complex  resonant  frequencies 
are  calculated  for  different  configurations  versus  substrate  height  h2.  The  convergence  of 
the  results  is  investigated  and  ensured  by  using  a  proper  number  of  basis  functions.  It 
is  shown  that  there  are  dual  resonant  frequencies  of  the  stacked  structure  for  each  mode 
which  are  related  to  the  constitutive  resonators  of  the  structure  and  a  function  of  the 
mutual  coupling  between  them.  Throughout  the  analysis,  the  exp(— iu>t)  time  dependence 
is  used  and  suppressed. 


2.  DYADIC  GREEN’S  FUNCTION  FORMULATION 

Consider  the  horizontally  stratified  medium  shown  in  the  Figure  2  where  each  layer  l  is 
of  thickness  hi  and  characterized  by  permittivity  ej  and  permeability  /ij.  The  two  coaxial, 
circular  conducting  disks  are  carrying  electric  current  distributions  Ji(r')  and  J j(r')  and 
located  in  medium  1  and  medium  2,  respectively.  Configurations  A  and  B  are  obtained, 
respectively,  by  placing  Ti(r')  at  =  do  and  z[  =  d\. 

The  electric  fields  in  regions  1  and  2  due  to  the  current  distributions  and  J^r*) 

may  be  expressed  as 

JffYW'ZiAW)-7i<?)  +  ™l*  J^dV'SliJ(r,F').  J2(rl)  (la) 
Z,(F)  =  u>px  JJfYt  dV'V^r,?)  •  7i(F/)  4-  JJJvi  dV'Uit2(r,r')  •  71(f')  (16) 

where  r1)  is  the  dyadic  Green’s  function  relating  the  electric  field  in  region  /  to  an 

impressed  current  source  in  region  p. 

For  the  stacked  microstrip  antenna  structure,  the  current  distributions  are  surface 
currents  at  z‘  =  z[  and  z'  =  z'2 


JiO1')  =  V)*(*'  -  z[) 


(2a) 
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H?)  =  Jf'VW*'  -  4) 


(26) 


Using  a  dyadic  Green’s  function  formulation  for  horizontally  stratified  media  [21]- 
[23],  we  may  write  the  spectral  domain  dyadic  Green’s  function  Gip(r,r*)  in  cylindrical 
coordinates  as 


I  m=  — OO  ® 


(3) 


where  $p\p,  z\p,  z')  is  given  in  the  following  for  the  four  cases  of  source  and  observer 
positions.  When  both  observation  and  source  positions  are  in  region  l  with  z  >  z' ,  we  have 

J&i 0. *;?'.*')  =  RT^Tp^  + 

+  j*V  l  -  R™ ^TMf, ku)e‘klMt  +  Jlffie  3  '*  1  'vm(/>,-ku)e  '•  ] 

k(,'. 

(4) 

where  fc/x  satisfies  =  fcf  —  A:*.  When  both  observation  and  source  positions  are  in  region 
/  with  z  <  z',  we  have 

if, ?'.*')  = 

[a™.®-*-  i'Jp',  k,.)rik‘-‘‘  +  h'jp', 

Me'*-'  + 

(5) 

When  the  observation  point  in  region  l  is  above  the  source  point  in  region  p,  we  have 

K;Wi?V)  =  [M/>, *..>«*■•*  + 


vTM 

AuJ 


fcpX  1  -  R™R™ei3h*‘K’ 


(6) 
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When  the  observation  point  in  region  l  is  above  the  source  point  in  region  p,  we  have 


1  X™  r  ■ 

+  1  -  R™RT [^1  e  **  ,s  +  vm(p,  — fc/«)c  tfc,,r] 


fcnp  •*cnp 


(7) 

In  the  previous  equations  (4)-(7),  the  dyadic  Green’s  function  is  decomposed  with  TE  and 
TM  wave  components  with  the  following  unit  vectors 


hm(P) 

h'miP  1  ±fc/x) 


pj^MKp)  -  4>J'm(k,P) j 


(8a) 

(86) 


6m(p»i^/i)  — 

vm(p\±ku)  = 


^p-j^J'miKp)  t  ^YxVLpJm{<kpp^ + 

TP^JUKp')  ±  4'j;  j^Mk*')  -  ^Jm(Kp') 


(9a) 

(96) 


•  It  may  be  noted  that  the  unit  vector  hm(p,  ±fcj*)  is  in  the  direction  of  the  electric  field  for 

a  horizontally  polarized  TE  wave  and  that  the  unit  vector  vm(p,  ±kXt)  is  in  the  direction 
of  the  electric  field  for  a  vertically  polarized  TM  wave. 

The  reflection  coefficients,  R^i  and  R%X,  at  the  upper  and  lower  boundaries,  respectively, 
^  of  layer  l,  and  are  given  by  the  following  recursion  relations 


„  flTn-i)  +  ay, 
n‘  l  + 


(70a) 

(106) 
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where  R£q  =  R£q  =  0  and  R*f  =  -1  and  =  1.  The  Fresnel  reflection  coefficients 
Rj(f±i)  and  are  defined  by 


pTE 

nl(l±\) 


pTM 

nt(l±l) 


f*(i±i)ku  -  mk(i±i)t 
P(l±i)ku  + 

t(l±l)ktx  -  Ctfyi l)x 
e(t±i)kix  +  ej^(i±i)* 


(11a) 

(11  b) 


The  transmission  coefficients,  X and  X™,  for  observation  positions  above  the  source 
points  are  given  by  the  following  recursion  relations 


e'V+o.4  [l  +  jffq+i)]  -E 
e'h<  d‘  [1  4-  R£fei2k*>h»]  Au(!+i) 

kt  k{l+l)teih<‘^d‘  [l-^ij  TM 
Jfc(/+1)  ku  eih<‘d*  [1  - u(,+1) 


(12a) 

(126) 


for  l  =  0, 1, ...  ,p  —  1  where  X™  =  X™  =  1  and  those  for  observation  positions  below  the 
source  points,  X™  and  X™,  are  given  by  the  following  recursion  relations 


■TE 

«-*"*  [i  +  «ff]  yTE 

(13a) 

n(J+ 1)  — 

•TM 

V 1)  *<*  [l  "  *3*]  ytm 

(136) 

n(i+l)  “ 

ki  k{l+ !),  _  /2TM+i)ei2fc(,+l)lfc(l+l)j  ni 

for  l  =  p  +  l,p  +  2, ...  ,3  where  X™  =  X™  =  1* 

The  transverse  (to  z)  components  of  the  spectral  domain  dyadic  Green’s  function 
may  be  written  in  the  following  form 

p.,0‘.'')]T  =  £;  £  C  Kp)  -?,,(*>„  z,z')  (14) 

m=  —  oo 

asf  t 

where  Jm(kpp)  is  the  complex  conjugate  transpose  of  the  kernel  of  the  vector  Hankel 
transform  (VHT)  [24]  given  by 


Jm(k„p) 


_  '  J'mM  i%Jm(kpp) 
.  irMKp)  j'm(kpP) 
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and 


zi  z  )  ~ 


'TM 


0  I 

L  o  CT/J 


For  source  and  observer  in  region  1  (assuming  z  >  z'),  we  have 


>TM 

‘.1,1 

»TE 


*? 


1  - 


Cl  C  _ 

1,1  — 


kiz 


1  - 


(16) 

(17a) 

g‘fcl  ,(*-*') 

(176) 

where  fcj,  satisfies  the  dispersion  relation  kft  +  k*  —  k2  =  For  an  observer  in  region 

2  and  source  in  region  1,  we  have 


/*TM  _ 
‘>J.i  — 


CTE  - 
^J.i  _ 


[i  -  *™j  [i  -  [i  -  sav”"1*-'1]  ,i  M_ 

k‘  [l-ll™fl™ei“'-k‘][l * 

1  [l  +  fl™]  [l  +  [l  +  fl™ 

k„  [1  -  [1  +  fl7fe‘“>'k>]  * 


(18a) 


(186) 


For  source  and  observer  in  region  2  (assuming  z  >  z'),  we  have 


/*TM  _ 
‘•J.J  ~ 


*2 


1  - 


-TE 


<’,-J  k u  1  -  R^fR^fei3k*^ 


-*)1 

(19a) 

<h)] 

J  c‘*»*(*— ') 

(196) 

For  an  observer  in  region  1  and  source  in  region  2,  we  have 


/*TM  _ 
>1.2  “ 


l  fl  -  ff™  fl  -  ptm  fl  -  pTMi2fc,.(.' -<*,)] 

ku  l1  fua  \[L  7tm  e _ \[l  ^ 

Jfef  [1  -  [!  _ 


/TE  _ 
>1,3  _ 


i  [i  +  ifff]  [i  +  «3f «»>■•<*->]  [i  -i  jg» ««».■('-*>] 


(20o) 


‘i. 


[1  -  I&RZf  [1  +  i£fe<l*»«*‘] 


c»*l  *  ( t-di )  gifc»,  (<<i  -  *' ) 


(206) 
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Due  to  reciprocity,  Pi(2il(kp,z,  z')  =  p2(^2(kpyz\z)  where  a  =  TE,TM. 

The  transverse  components  of  the  electric  fields  for  the  mth  mode  at  z  =  z[  and  z  =  z'2 
vanish  on  the  conducting  disks  and  may  be  expressed  as 


[Ei>m(p,z  =  z[ 

)1  = 

'Jr  2 

J  dkpkpJ m(kpp)  •  Ci  i  •  #!,*(&„) 

ujp2 

2 

JQ  dkpkpJ m(kpp)  -Ci(2  x<=x,  4=t, 

=  0, 

P  <CL\ 

(21a) 

[^l,m(/>i  Z  =  22j 

>1  =-^ 
Jr  2 

f  dkpkpJm(kpp)  •Cj.i)  ,  ,  .  ' 

“Pi 

2  , 

JQ  dkpkpJ m ( kpp)  -C2)2  x,=x^x_^ 

=  0, 

P  <a  i 

(216) 

where  TC^(kp)  ( j  ==  1,2)  is  defined  by 

*!?(*.)  =  jT  dppfjk.p)  ■  K^(p)  (22) 

and  where  K^(p)  is  the  Fourier  coefficient 

T?%)  =  ^  J*  (23) 

With  the  use  of  the  vector  Hankel  transform,  we  obtain 

TQ\p)  =  f~dpp7m(k,p)-TQ\kt)  =  0,  p>a ,  (24a) 

TQ\p)  =  /o“  dppl„(k,p)  ■  T?’](k,)  =  0,  p>  a,  (244) 

Equations  (21)  and  (24)  constitute  a  set  of  coupled  vector  integral  equations  for  the 
currents  K^\kp)  and  K^\kp).  For  configuration  A,  z{  =  dx  =  -hi  and  z'2  =  <f2  = 
-(hi  +  h2),  and  for  configuration  B,  z[  =  do  =  0  and  z2  =  d2  =  -(hi  4-  h2). 


3.  GALERKIN’S  METHOD 


Galerkin’s  method  is  now  employed  to  solve  the  coupled  vector  integral  equations  of 
q  (21)  and  (24).  The  currents  on  the  circular  disks  are  expanded  in  terms  of  the  complete 

orthogonal  set  of  TE  and  TM  modes  of  the  corresponding  cylindrical  cavities  with  magnetic 
side  walls 


where  j  =  1,2  and 


=  E  +  E 

n=l  p=l 

*!?(/>)  =  E  «!25S(p) + E  Ki&lw 


r  =  l 


»=1 


?!£(/>)  = 


JmiflmnP/Q’j) 

0, 


,  for  p  <  aj 
for  p  >  aj 


(25  a) 
(25  6) 


(26a) 


^mp(p)  = 


^Jm(ampPl  &j) 


»*>/> 

J!n{ampp/<lj) 

0, 


(266) 


,  for  p  <  a j 
for  p  >  aj 

flmn  and  amp  are  the  nth  and  pth  zeros  of  J„(0mn)  =  0  and  Jm(amp)  =  0,  respectively.  In 
the  above,  N  and  P  correspond  to  the  number  of  TM  and  TE  basis  functions,  respectively, 
taken  for  the  upper  disk  and  R  and  5  correspond  to  those  taken  for  the  lower  disk. 

The  corresponding  vector  Hankel  transform  of  the  currents  K^(kp)  ( j  =  1,2)  are  given 
by  the  following 


where 


*£’(*,)  =  E  +  E 

n=l  p= 1 

*!?(»,)  -  E  +  E 


r  =  I 


»=1 


££(*.)  =  An 


-  k* 


(27  a) 
(276) 


(28a) 
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T<j)/i  \  _  J'm  ( )  0 

^  kl-ia^/atf  [Jm(*pa;). 

Substituting  these  expressions  into  equation  (21),  we  obtain  for  the  mth  mode 

=  2l)]r  =  E  amn  /  dkpkp  Jm( fcpp)  •  j  *  V’mn(^p) 

1  n=1  *  0  *  *i 

+  E‘S  -?2J(fc.) 

p=1  •'0  ix-Zj.x-Zj 

12 

+  E«S  /"  dk,k,7„(k,?)  ■  -iS(M 

r=1  •'O  ix-Xj.z-Zj 

5 

■b  E  ^m*  J  dkpkpj m(kpP)  •  ,2 1 

,=l  Jo  Iz-Zj.z-z, 

=  0,  p  <  aj 

[E,.m(?,  z  =  4)1  =  E  «SJi  r  dk,k,7m(k,p)  ■  ?,tl|  .  • 

Ji  n=1  *0  i*  — *&i«— 

p 

+  E  “4  r  dk,k,7„(k,p)  ■  ?mL  .  .  •  ?£>(*.) 

p=1  ^0 

/l 

■b  XI  amr  /  dkpkpj m(kpP)  •  (j,a| 

+ e  r  dkpkpim(kpP) .  j2ti  |  .  $*\(kp) 

=  0,  p  <  aa 


(286) 


(29a) 


(296) 


where  we  have  defined  Z[  p(kp}z,z')  =  —  ^Ci^ikp,  z,  z'). 

Multiplying  (29a)  by  pV’iii*  (p)  &nd  p^j\p)  and  integrating  from  zero  to  infinity  for 
h  =  1, 2, . . . ,  N  and  *  =  1, 2, . . . ,  P,  and  applying  Parseval’s  relation  for  the  vector  Hankel 
transform  [24],  we  obtain 

E  °mn  /  dkpkpj (kp)  •  |u|  ’  ^mn(K) 

+ E  >4  f  dk,kj£"(k,)  •  ....  ■  ?!£(*,) 

P_J  <0  ix-Zj.z— z, 

il 

+  E  °mr  /  dkpkp^(kp)  •  ?u|  • 

+ E  «£J  /“  dk,k,^!(kt)  ■  .  •  ?£(*,) 

4=1  Jo  '*  -**•*-*! 


h  =  1,2,. ..,1V 


(30a) 
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5Z  °mn  f  dkpkpfimi  (kp)  '  fl.l  ,  '  'i’Llii^p) 

+  51  &mp  J  dkPkP^Ji  (K)  1  Iwl^  ,  •  ^ttp(fc/>) 

P=1  ,#0  —  X|.X“  *1 

+  T,amr  f  dkpkp$^(kp)  •  f,t,|  •  ^J(A:p) 

+  5Z  J  dkf>kP<t>mi  ikp)  ‘  fl.i)  ,_  ,  ,  *  fimlikp) 

*=  1  '*-*>»*-*! 


t  =  1,2, ...,P 


(306) 


Similarly,  multiplying  (29b)  by  and  p0^^(p)  and  integrating  from  zero  to  infinity 

for  j  =  1,2,...,/?  and  k  =  1 , 2, . . . ,  5,  and  applying  Parseval’s  relation,  we  obtain 

e  «si  r  dkfk^(k,)  ■  ?M|  •  #i(k,) 

n_i  •'0  IX-*,,*— 

+  5Z  ^mp  JQ  dkPkP*l>mj  ikp)  ‘  fj,l  x,_  ,  x  ,  *  $mp(kp) 

P—  1  » 1  *3 

+  H  aLr  J  dkPkp$m]  (kp)  *  ,  .  .  *  ^mr(^p) 

+ £  e  r  -MkCtw  •  ,_4  •  «2ow 

«=1  ,/0  '*  -*»'*-*» 


j  —  1,2,...,/? 


(30c) 


°  =  S>mi  f  <‘KK^(K)  •  {,,  |  ■  ?£>(*,) 

+  Hbmp  f  dkPkp$mk  (M  -  .  *  ^mp{kp) 

+  51  amr  J  dkpkpfirX  (kp)  *  (a,a|  ,  .  *  V’i.r^p) 

+  5Z  Jo  dkPkp$mk(kp)  •  •  ^ml{kp) 


j  —  1 » 2 , . . . ,  S 


(30  d) 


These  four  equations  (30a)-(30d)  constitute  a  system  o{N  +  P  +  R  +  S  linear  algebraic 
equations  which  may  be  written  in  matrix  form 

1-2  =  0 
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(31) 


where 


'mCT 

L^Ijvxp 

r 

Ia\NxR 

L^Jpxp 

l-^JPxP 

[A]*"*" 

Ia1pxs 

l^J/lxAT 

L-^JPxp 

I-^JPxP 

fAl***^** 

1-aJPx5 

l>1J5xJV 

r  41^2>^1) 
l^JSxP 

K*JSxP 

r4i^aM,) 

K*JSxS 

rteiv.i 

MJW 

l  M2>]ix.  J 

Each  element  of  the  submatrices  of  A  is  given  by 


AfkX<,)  =  fQ  ^KKirrlliK)  • 


(32) 


(33) 


(34) 


7$n(kp)  J$n(kp)  represent  either  i^n{kp)  or  4>£n(kp).  More  explicitly,  the 
integrals  of  (34)  may  be  written  as 


_  _  j  \  a  j  (a  \  /"  jl  l  _ ^m{kPal)'^mi.kPav) _ fTM/j,  \ 

'*  ‘  J"  "  [(/W*)’  -  *;]  [(au/M*  -  *j]  "  ' 


(35a) 


A»  -~Ati  -xma,ar~0~' 7">(a"‘)/0  *«*'  -  *j  {"(*') 

’  ”  ’  [(ami/a,)>  -  *j]  [(<W«,)'  -  *:j]  ,J> 


As  the  resonant  frequencies  are  complex  due  to  radiation  loss,  the  branch  point  and 
the  surface  wave  pole  singularities  of  the  functions  p )  (a  =  TE,  TM)  are  also  complex. 

Therefore,  the  path  of  integration  in  the  complex  plane  is  deformed  below  the  real  axis 
to  avoid  these  singularities.  The  integrals  of  equations  (35)  are  evaluated  numerically  along 
the  integration  path  shown  in  Figure  3. 


For  nontrivial  solutions  of  the  system  of  equations  (31),  it  is  required  that 


det  |A(w)|  =  0 


(36) 


This  is  the  eigenvalue  equation  for  the  stacked  circular  microstrip  antenna.  Since  the  basis 
functions  (25)  form  a  complete  set,  the  eigenvalues  of  (36)  correspond  to  the  exact  resonant 
frequencies  of  the  stacked  microstrip  antenna  structure  as  N ,  P,  R ,  and  S  — ►  oo.  With  a 
judicious  choice  of  basis  functions,  usually  only  several  basis  functions  are  needed  to  obtain 
accurate  results. 


4.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  complex  resonant  frequencies  in  the  following  results  are  calculated  versus  the 
parameter  /ij/a2,  which  from  a  practical  viewpoint,  is  the  parameter  most  easily  varied  if 
the  substrates  of  medium  1  and  medium  3,  on  which  the  microstrip  disks  are  printed,  are 
fixed.  Attention  is  given  to  the  lowest  order  mode  with  m  =  1.  This  mode  is  the  lowest 
order  mode  which  radiates  maximally  in  the  broadside  direction  and  is  widely  used  in 
microstrip  antenna  applications.  For  the  lowest  m  —  1  mode,  convergence  of  the  solution 
for  the  different  configurations  presented  is  confirmed  by  varying  the  number  of  current 
basis  functions  on  the  upper  and  lower  disks,  which  are  taken  to  be  [(IV  =  2,  P  =  0),  ( R  — 
1,5  =  0)],  [(JV  =  2, P  =  1  ),(*  =  2,5  =  1)],  [(N  =  3,P  =  2),(R  =  2,5  =  2)],  and 
[(IV  =  4,P  =  3),(P  =  3, 5  =  2)j.  It  is  found  that  the  solutions  converge  and  that  the 
resonance  curves  using  basis  functions  with  [(IV  =  4,  P  =  3),  (R  =  3, 5  =  2)]  follow  closely 
those  using  basis  functions  with  [(jV  =  3,P  =  2),(R  —  2, 5  =  2)]  for  the  cases  presented, 
where  the  maximum  relative  difference  between  the  two  curves  is  less  than  one  percent. 
The  convergence  of  the  solution  for  the  lowest  m  =  1  mode  is  illustrated  in  Figure  4  for 
configuration  A  with  oi /aj  =  1,  hi  =  =  O.laj,  =  2.65cot  and  €j  =  l.Oleo* 

From  Figure  4,  it  is  shown  that  the  stacked  circular  microstrip  structure  has  dual  resonant 
frequencies  associated  with  the  two  disks.  As  illustrated  in  the  figure,  the  number  of  basis 
functions  must  be  sufficient  so  that  the  correct  resonant  behavior  near  hj/aj  =  0.8  is 
obtained,  or,  so  that  the  coupling  interaction  between  the  two  disks  may  be  adequately 
accounted  for. 
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To  explain  physically  the  resonant  behavior  of  the  stacked  microstrip  structure,  it  will 
be  shown  that  the  dual  resonant  frequencies  may  be  related  to  the  resonant  frequencies  of 
two  constitutive  resonators  of  the  structure  which  axe  electromagnetically  coupled  through 
their  fringing  fields.  For  thin  substrates,  the  constitutive  resonators  have  resonant  frequen¬ 
cies  close  to  those  of  cavities  with  electric  top  and  bottom  walls  and  lossy  magnetic  side 
walls.  Considering  configuration  A,  one  resonance  of  the  stacked  structure  is  associated 
with  the  cavity  formed  by  the  lower  disk  and  the  ground  plane  illustrated  in  Figure  5(a). 
The  second  resonance  is  considered  in  two  limits  of  the  substrate  height  h2.  In  the  small 
hi  limit  when  the  ratio  aj/aj  is  not  too  large,  the  second  resonance  may  be  associated  with 
the  cavity  formed  by  the  two  disks  as  shown  in  Figure  5(b).  For  large  h2,  the  second  res¬ 
onance  may  be  associated  with  the  cavity  formed  by  the  upper  disk  and  the  ground  plane 
as  shown  in  Figure  5(c).  In  comparing  the  resonances  of  the  stacked  microstrip  structure 
to  those  associated  with  the  constitutive  cavity  resonators,  the  resonances  associated  with 
the  cavities  of  Figures  5(a)  and  5(c)  are  calculated  using  the  configurations  of  Figures  6(a) 
and  6(c),  respectively.  The  resonant  frequency  associated  with  the  cavity  of  Figure  5(b) 
formed  by  the  two  disks  is  calculated  approximately  using  the  configuration  of  Figure  6(b). 
Configuration  B  may  be  considered  in  a  similar  manner. 

To  illustrate  the  relation  between  the  resonances  of  the  stacked  microstrip  structure 
to  those  of  its  constitutive  resonators,  shown  in  Figure  7  are  the  resonant  frequencies  for 
the  lowest  m  =  1  mode  of  configuration  A,  compared  with  the  resonances  associated  with 
the  configurations  of  Figure  6.  The  open  circles  and  squares  on  the  solid  curves  are  used 
to  help  distinguish  the  two  pairs  of  resonant  frequencies.  For  small  fi2,  the  two  resonances 
correspond  to  ksa 2  =£  and  fc2a i  =  flu  where  /?u  =  1.84.  It  is  seen  that,  for  comparable 
disk  radii,  the  ratio  alA//i2«2/a2v//i3C3  approximately  determines  the  separation  of  the  two 
resonant  frequencies  for  small  hi.  It  is  shown  that  in  both  the  small  and  large  hi  limits,  the 
resonant  frequencies  approach  those  associated  with  the  constitutive  resonators.  Thus,  the 
coupling  interaction  between  the  two  constitutive  resonators  in  these  two  limits  is  small. 
For  intermediate  values  of  /i2,  the  resonant  frequencies  associated  with  the  constitutive 
resonators  are  affected  by  the  electromagnetic  coupling  between  them.  As  in  Figure  7, 
the  two  resonance  curves  essentially  follow  those  of  the  constitutive  resonators  with  the 
most  deviation  occuring  where  the  real  parts  of  the  resonant  frequencies  of  the  constitutive 
resonators  are  approximately  equal.  The  presence  of  the  upper  disk  increases  that  radiation 
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loss  associated  with  the  cavity  formed  by  the  lower  disk  and  the  ground  plane. 

To  illustrate  the  effect  of  different  disk  radii  on  the  resonance  frequencies,  shown  in 
Figures  8  and  9  are  the  resonant  frequencies  for  configuration  A  with  ax/a2  =  0.9  and 
°i/oj  =  1.1,  respectively.  The  constitutive  resonators  with  aj/oj  =  0.9  are  less  coupled 
than  those  with  ax/a2  =  1.  This  is  consistent  with  the  resonance  curves  in  that  the  resonant 
frequencies  with  ax/a2  =  0.9  follow  more  closely  the  constitutive  resonances  than  those  of 
Figure  7  where  ax/a2  =  1.  Conversely,  the  resonators  with  ax/a2  =  1.1  are  more  tightly 
coupled  than  those  with  ax/a2  =  1,  and  therefore,  the  resonance  curves  deviate  significantly 
from  the  constitutive  resonances  in  the  region  of  maximum  coupling  interaction.  In  fact, 
in  the  case  ai/a2  =  1.1,  the  coupling  is  such  that  in  the  region  where  the  constitutive 
resonances  are  approximately  equal,  the  resonances  of  the  coupled  structure  split  such 
that  the  resonant  curves  follow  opposite  constitutive  resonances  in  the  small  and  large  /i2 
limits. 

The  effect  of  increasing  the  substrate  heights  is  illustrated  in  Figure  10  where  the 
resonance  curves  for  configuration  A  with  ai/o2  =  1  and  hi  —  hs  =  0.2a2  are  plotted.  The 
coupling  between  the  resonators  is  tighter  in  this  case  relative  to  that  of  Figure  7  where 
=  0.1o2  due  to  the  increase  in  fringing  fields. 

The  resonant  frequencies  for  configuration  B  with  ai/a2  =  1,  hi  =  hs  =  0.1a2,  cx  = 
e3  =  2.65eo,  and  e2  =  1.01eo  are  shown  in  Figure  11.  In  this  case,  as  h2  approaches  zero, 
both  resonant  frequencies  approach  the  same  value.  However,  due  to  the  presence  of  the 
air  gap,  when  h2  is  small,  the  resonant  frequency  associated  with  the  cavity  fromed  by 
the  two  disks  is  a  rapidly  varying  function  of  h2,  making  this  configuration  undesirable  in 
some  practical  cases  where  height  tolerances  cannot  be  maintained. 

The  advantage  of  the  stacked  structure  over  conventional  microstrip  antennas  is  the 
introduction  of  the  second  resonance  which,  when  close  to  the  first  resonance  may  provide 
for  wide  bandwidth  operation,  or,  when  separated  from  the  first  resonance,  may  provide 
for  dual  frequency  operation.  In  addition,  the  radiation  loss  associated  with  the  cavity 
formed  by  lower  disk  and  the  ground  plane  may  be  enhanced  by  the  presence  of  the  upper 
disk,  making  the  stacked  microstrip  antenna  a  better  radiator  than  conventional  microstrip 
antennas. 


5.  CONCLUSION 


The  complex  resonant  frequencies  of  the  open  structure  of  a  microstrip  antenna  con¬ 
sisting  of  two  circular  microstrip  disks  in  a  stacked  configuration  have  been  rigorously 
calculated  as  a  function  of  the  layered  parameters  and  the  ratio  of  the  radii  of  the  two 
disks.  Using  a  dyadic  Green’s  function  formulation  for  horizontally  stratified  media  and 
the  vector  Hankel  transform,  the  mixed  boundary  value  problem  is  reduced  to  a  set  of 
coupled  vector  integral  equations.  Employing  Galerkin’s  method  in  the  spectral  domain, 
the  complex  resonant  frequencies  are  calculated  and  convergence  of  the  results  is  demon¬ 
strated.  It  is  shown  that  for  each  mode,  the  stacked  circular  microstrip  structure  has  dual 
resonant  frequencies  which  are  associated  with  the  two  coupled  constitutive  resonators  of 
the  structure  and  which  are  a  function  of  the  mutual  coupling  between  them.  This  mutad 
coupling  depends  on  the  geometrical  configuration  of  the  stacked  structure,  the  layered 
parameters,  and  the  disk  radii.  In  the  large  and  small  h 2  limits,  the  resonant  frequencies 
of  the  stacked  microstrip  structure  approach  those  of  the  constitutive  resonators.  The 
maximum  coupling  effect  occurs  where  the  real  parts  of  the  resonant  frequencies  of  the 
constitutive  resonators  are  approximately  equal,  where  the  behavior  of  the  resonances  in 
this  region  is  a  function  of  the  coupling.  The  dual  frequency  behavior  of  the  stacked  mi¬ 
crostrip  structure,  easily  controlled  by  varying  the  parameters  of  medium  2  and  disk  radii 
ratio,  given  fixed  parameters  for  medium  1  and  medium  3,  may  be  used  to  broaden  the 
bandwidth  or  provide  for  dual  frequency  use  of  the  antenna. 
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PROPAGATION  PROPERTIES  OF  STRIPLINES 
PERIODICALLY  LOADED  WITH  CROSSING  STRIPS 
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Department  of  Electrical  Engineering  and  Computer  Science 
Massachusetts  Institute  of  Technology  Cambridge,  MA02139 

Abstract 

A  rigorous  dyadic  Green’s  function  formulation  in  the  spectral  domain  is  used  to  study 
the  dispersion  characteristics  of  signal  striplines  in  the  presence  of  metallic  crossing  strips. 

A  set  of  coupled  vector  integral  equations  for  the  current  distribution  on  the  conductors 
is  derived.  Galerkin’s  method  is  then  applied  to  derive  the  matrix  eigenvalue  equation  for 
the  propagation  constant.  The  dispersion  properties  of  the  signal  lines  are  studied  for  both 
cases  of  finite  and  infinite  length  crossing  strips. 

The  effects  of  the  structure  dimensions  on  the  passband  and  stopband  characteristics 
are  investigated.  For  crossing  strips  of  finite  length,  the  stopband  is  mainly  affected  by  the 
period,  the  crossing  strip  length,  and  the  separation  between  the  signal  and  the  crossing 
strips.  For  crossing  strips  of  infinite  length  carrying  travelling  waves,  attenuation  along 
the  signal  line  exists  over  the  whole  frequency  range  of  operation. 
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I. Introduction 


In  microelectronic  computer  packaging,  a  problem  of  practical  interest  is  the  study 
of  propagation  characteristics  of  microstrip  lines  embedded  in  a  layered  medium  in  the 
presence  of  periodic  crossing  metallic  strips. 

The  analysis  of  striplines  and  finlines  with  periodic  stubs  has  been  studied  by  Kitazawa 
and  Mittra[l],  where  a  technique  based  on  the  network- analytical  formulation  is  used.  A 
slow-wave  coplanar  waveguide  on  periodically  doped  semiconductor  substrate  has  been 
carried  out  by  Fukuoka  and  Itohi2l.  Gu  and  Kong[3]  used  a  quasi-static  approach  to  study 
single  and  coupled  lines  with  capacitively  loaded  junctions.  The  propagation  characteristics 
of  signal  lines  in  a  mesh-plane  environment  has  been  presented  by  Rubin[4].  More  recently, 
the  propagation  characteristics  of  signal  lines  in  the  presence  of  periodically  perforated 
ground  plane  is  studied  by  Chan  and  Mittra[5j. 

An  analysis  of  a  width-modulated  microstrip  periodic  structure  using  quasi-static  ap¬ 
proach  is  presented  in  [6].  A  hybrid  spectral  domain  analysis  for  similar  periodic  structures 
have  been  carried  out  in  [7], 

In  this  paper,  hybrid  mode  analysis  is  used  to  study  the  propagation  characteristics  of 
striplines  periodically  loaded  with  crossing  metallic  strips.  The  periodic  crossing  strips  are 
assumed  to  have  finite  or  infinite  length.  A  dyadic  Green’s  function  formulation  for  the 
periodically  loaded  structure  is  derived.  A  coupled  set  of  vector  integral  equations  for  the 
surface  current  distribution  is  formulated.  Galerkin’s  method  is  then  applied  to  transform 
the  resulting  set  of  integral  equations  for  the  current  distribution  into  a  determinantal 
equation  from  which  the  dispersion  characteristics  are  obtained. 

The  propagation  properties  of  one  signal  line  and  two  coupled  lines  in  the  presence  of 
periodic  crossing  strips  are  investigated.  Numerical  results  for  the  passband  and  stopband 
characteristics  are  presented. 


2 


II. Dyadic  Green’s  Function  Formulation 


The  geometrical  configuration  of  the  problem  is  shown  in  Fig.l  where  M  signal  striplines 
located  at  z  =  zm ,  m  =  1, 2,  •  •  • ,  Af,  are  periodically  loaded  with  crossing  metallic  strips 
having  a  period  p.  The  crossing  strips  are  of  width  wc,  length  Lc ,  and  axe  located  in  the 
plane  z  =  xa^+i.  Both  the  signal  lines  and  the  crossing  strips  are  embedded  in  the  same 
layer  (!)  having  parameters  (c(,/i0)- 


In  general,  the  electric  field  can  be  expressed  in  terms  of  the  dyadic  Green’s  function 
and  the  current  distribution  on  the  strip  surfaces  as[8] 

Ei(r)  =  luip,  JJJ^  Gu(fy)  ■  J{f,)dV'  (1) 

where  Gu(f,  f*)  is  the  dyadic  Green’s  function  when  both  the  observation  point  r  and  the 
source  point  f1  are  located  in  the  Ith  layer  of  the  stratified  medium.  For  z  >  z' , 

(1  -  €>«*•>*)  +  -Rif  ( 

U(klt)e-ih •'*«’  +  R^v(-ku)etkt‘z> }  ). 


( 


For  z  <  z', 


{ (T  . h-  sff 

[h(-klz)e'hl‘<  +  Rif  A(jfeu)e<fc,*(a*,_i;)] 

+  (1  -  +R^v{ku )etk'- *'] 

v{-ku)elh,‘t<  +  Cp(lt|J)e’*,-(n,-,!)]  j 


(2a) 


(26) 
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where  n  and  zj  axe  the  local  coordinate  defined  as  zj  =  z  4-  di,  z\  =  z'  4-  dj,  and 
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kg  —  ikz  4"  yky 


kj  =  k]  4-  k]z,  kt  =  \kt 


h(±kit) 

v(±kiz) 

r. 


zky  y  k  2 
kizkg  *  kg 

*~kiK  +Z~ki 

iz  4- yy,  rt-xx+yy. 


(3) 


In  (2),  iZjf*  and  R^f  are  the  reflection  coefficients  of  TM  mode  and  TE  mode  at  the 
upper  boundary  of  the  /th  layer,  and  R^f  are  the  reflection  coefficients  of  TM  mode 
and  TE  mode  at  the  lower  boundary  of  the  ith  layer.  They  cam  be  obtained  recursively  as 


RZi  = 


l  +  RT(i- 


a  =  (  TE,  TM  ) 
a  =  (  TE,  TM  ) 


(4.) 

(46) 


where  and  axe  the  Fresnel  reflection  coefficients  of  a  mode  across  the  inter 

faces  at  z  =  — of|— i  and  z  =  —  dj,  respectively.  The  explicit  forms  are 


nTS  _  ^ lx  fy±i)» 

ku  4-  k(i±i),  ’  eiiifcj*  4- c/A:(i±i), 


For  our  problem,  only  transverse  currents  J,  having  no  z  component  exist,  the  transverse 
electric  field  (to  z)  Ei ,  in  layer  (/)  is  thus  given  by 


Eig{r)  =  iu>n j  JJ^drtGji(f,rg)  ■  Jt{r,)dS' 


(6) 


where  Gjl^f,?)  is  the  (2  x  2)  transverse  (to  z)  part  of  the  dyadic  Green’s  function,  and 
can  be  expressed  in  the  k,  domain  as 


m 
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where  for  z  >  z'} 


The  transverse  electric  field  can  be  expressed  using  Floquet  harmonic  representation  in 
the  y  direction  as 

r°°  fP/2 

Ei,{r)  =  ium  dx'  dy'Gflp{r  (9) 

J-oo  J-p/2 


where  is  given  by 


where  Jfeyn  =  Jfey0  -I-  2nir/p,  ky0  is  the  propagation  constant  of  the  dominant  harmonic  in 
the  Floquet  representation.  We  assume  that  we  have  M  signal  striplines  and  one  crossing 
strip  within  one  period.  Thus,  Jt  cam  be  expressed  as 

J/ )i  U>m/2  X  ^  Xm  ■+■  Wm/2 ,Z  =  Zm 

(H) 

0,  elsewhere. 
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Substituting  (11)  into  (9),  we  have 


M  +1  oo 


Etl(r)  =  tuw(2w)3  £  dk*eihm*fii(kz,kyn,z,zm)  •  jfeyn)  (12) 

m=l  n=  — oo  *  —co 


where 


rp/J 


Jm(K,k yn)  =  ~  [  dx'e  **•*'  f  dy'e  ih’nV' Jm(z' ,y'). 

Zlr P  ^ -oo 


(13) 


The  electric  field  Eu  given  by  (12)  satisfies  the  boundary  conditions  at  the  interfaces 
between  the  dielectric  layers  of  the  stratified  medium.  Imposing  the  final  boundary  condi¬ 
tion  that  the  tangential  electric  fields  vanish  on  the  metallic  surfaces  of  the  signal  striplines 
and  the  crossing  strips,  we  get  a  set  of  vector  integral  equations  for  the  current  distribution 
on  all  the  metallic  strips.  Thus,  we  have 


Af +1  oo 


m=l  n=  — oo  J -oo 

*q  ~  w,/2  <  z  <  *,  +  w,/2,  -p/2  <  y  <  p/2,  z  =  z„ 

?  =  !,•••,  M  (14a) 


M+l  oo 


m=l  n=  — oo  *  00 


-  Tc/2  <  Z  <  Ic/2,  -u>e/2  <  y  <  we/2,z  =  zm+i  (146) 

where  (14a)  satisfies  the  boundary  condition  on  the  M  signal  strips,  and  (14b)  satisfies  the 
boundary  condition  on  the  crossing  strip.  The  task  is  to  solve  this  set  of  vector  integral 
equations  using  the  moment  method  to  get  the  dispersion  relation. 
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Ill-Numerical  Solution  for  the  Dispersion  Relation 

A)Qne  Signal  Strioline  Loaded  with  Crossing  Metallic  Strips 

In  this  section,  we  study  the  case  of  one  signal  stripline  in  the  presence  of  periodical 
crossing  strips  as  shown  in  Fig. 2.  The  signal  stripline  and  the  crossing  strips  are  located 
in  a  dielectric  layer  bounded  by  two  ground  planes  and  placed,  respectively,  at  z  =  z\  and 
z  =  Zj. 


To  apply  the  moment  method,  we  choose  an  appropriate  set  of  basis  functions  to  rep¬ 
resent  the  surface  current  Ji(z,y )  and  J2(x,y)  as 

.V,  JVj 

(x,y)  =  x  aiftj(x>y)  +  y  ^2  (15a) 

i=-S  t  ;=-ATi 

Nt 

h  (*,y)  =  i^2cjf3j(z,y)  +  y^dj/4j(*,y)  (156) 

;=i  j—0 

where  3\{zyy)  is  the  surface  current  on  the  signal  stripline,  J2(z,t/)  is  the  surface  current 
on  the  crossing  strip,  aj,  6;,  Cj,  and  dj  are  the  expansion  coefficients,  /i>(*,y),  /a j(*,y), 
f3j(x,  y)  and  /*j(*,y)  are  the  basis  functions.  The  explicit  forms  of  the  basis  functions  are 
as  follows  : 


where 


/a;(*,y)=T0(x,u;1)e<fc-» 
f3j(z,y)  =  Pj(z,Le)(p/2ir)T0(y,we) 
fa(x>y)  =  Tj(z,  Le)Q(y,  we) 


PM,v) 

m,T,) 


f  (l/*7)sin(2j'7r£/T7),  -rj/2  <t<i)/2 

0,  elsewhere 

f  (?/ cos(*{/tj),  —f)l2  <(<v/ 2 

0,  elsewhere 

f  cos(2 j  it  £/r))/  yJ{i)/2)2  -  {3,  -t?/2  <i<i)/2 


0, 


elsewhere. 


(16a) 

(166) 

(16c) 

(16d) 

(17a) 

(176) 

(17c) 
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The  surface  current  on  the  signal  stripline  is  basically  of  travelling  wave  type.  Due  to  the 
periodic  loading,  the  basis  functions  on  the  signal  stripline  are  chosen  as  a  superposition  of 
space  harmonic  modes.  On  the  crossing  strips,  the  surface  current  is  basically  of  standing 
wave  type,  and  the  phase  variation  along  the  y  direction  on  the  crossing  strips  can  be 
neglected. 


Let  Pj(kt,ri ),  (?(&„„,  77),  and  fj(kx,Tj)  be  the  Fourier  transforms  of  Pj(x,7}),  Q(y,Tj), 
and  Tj(x,rj)  respectively,  we  have 


=  (18o) 

<?(*,»,  1)  =  -  r  dyt-'i’-Q(y,V)  =  -  =  Q(-k,„,T,)  (186) 

P  J-fi/2  [kyn7])  T 


"7/2 

Tj(kz,Ti)  =  —  /  dxe-'kmZT)(x,T))  =  -[J0(kxr]/2  ~  jir)  +  J0(kxrj/2  +  jir)\  =  Tj(-kutrf) 

J- n/2  4 

(18c) 


where  Jo(a)  is  the  Bessel  function  of  the  zeroth  order,  Pj(kx,rj)  is  an  odd  function  of  kx , 
fj(kx,rj)  is  an  even  function  of  Jfe*,  Q{kyn,Tj)  is  an  even  function  of  kyn.  When  kxrj/2 
approaches  ±jir,  Pj(kx,r])  approaches  ±l/(47rt);  when  kynT]  approaches  ±tt,  Q{kvn,rj) 
approaches  l/(47r). 


With  these  basis  functions,  the  Fourier  transform  of  the  surface  current  Jm(kx ,  A:vn)  can 
be  derived  as 

N,  Ni 

J\{km,ky  n)  =  *  ^  1  ajPlj(kx,  kyn)  V  ^  1  8jF2j{kx,kyn)  (19a) 

j=-Nx  j=-Ni 

;=i  j=o 

where  Fij(kx,kvn)  is  the  Fourier  transform  of  /»>(*, y),  *  =  1,2, 3, 4.  The  explicit  forms 
are 


■F’lj(kX)  kyn)  —  8jnPl(kX,W\) 


(20  a) 
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F2  ]{kz,kyn)  —  o(^*>^l) 

F3j(kXfkyn)  =  fn(kyn,We}Pj(kzt  Le) 

F*j(kz,kyn)  =  Q(kyn,Wc)Tj(kx,  Lc) 


(206) 

(20c) 

(20d) 


where  is  the  Kronecker  delta  function.  Substituting  (20)  into  (14),  we  have 


X)  /  dkxeih'*gfl(kx,kyn,zuzl)-  [x  ^  ajFu(kmt  kyn)  +  y  bjF2j(kx,  fc„n)] 

n=  —  oo  J—oo  —  j—  —  Ni 

N, 

■t"Pj|(^*5  ^yni  )  -Zj)  •  ^  ^  ']cjFjj(kx,kyn)  ~t~  y  ^  >  djF^j{kx ,  fcyn)J  =  0, 


—  u>i/2  <  z  <  tci/2,  —p/2  <  y  <  p/2,  z  =  zi 


(21a) 


£  /°°«,ei‘-'^(A:„jb„„,2a,z1).[i  £  aiFli(k„k,n)  +  y  £ 

n=  —  oo  J~oo  j=— N\  j=—N\ 

r  n,  iv4 

+9u{kz,kyn,Z2iZ2)  •  jz  ^  ]  cjFzj(kx,kyn)  +  y  ^  ^  djFjj(kx,  fcyn)j  =  0, 


—  Lc/2  <  x  <  Ld 2,— ttfc/2  <  y  <  wc/2,z  =  z2. 


(2U) 


Applying  Galerkin’s  method,  we  choose  xe-,k*k*Pi(z,Wi)  (&  =  — JVj,--*, N2)  and 
ye-ih*k*To(z, w2)  (k  =  -Ni,--’,N2)  as  testing  functions  for  the  signal  stripline.  Tak¬ 
ing  the  inner  product  with  (21a),  we  obtain  2(N\  +  N2  +  1)  equations.  Similarly,  taking 
the  inner  product  of  £/ji(*>y)  (»  =  l,***,iVj)  and  y/4*(z,y)  {k  =  0, •••,JV4)  with  (21b), 
we  obtain  another  N2  +  N4  +  1  equations. 


After  arrangement,  we  get  the  following  matrix  equation  : 


r  7o.i) 

zij 

7(1.2) 

7(1.3) 

z«,; 

7(1.4) " 
ZiJ 

r°; 

7(2, 1) 
Zi,j 

7(3.2) 

zi.i 

7(2.3) 

7(2.4) 

Z<,; 

h 

7(34) 

7(3,2) 

7(3,3) 

7(3,4) 

z*.i 

ci 

7(4,1) 

L 

7(4,2) 

zi,i 

7(4,3) 

z<,; 

7(4,4) 

z»,i  J 

(22) 
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(23) 


Each  entity  Z\ rf)  in  (22)  is  a  submatrix,  the  explicit  form  of  the  elements  is 


E  dk,s\rf\k.,k 

_ _  J  —  oo 

n=  —  oo 


where 


SiJ9\ka,kyn)  —  f'ri(—kx,—kyn)gil^(hX)kyntZttzm)Fqj(kz,kyn) 

\zu  for  r  =  1,2  (  zit  for?  =  1,2 

=  <  ,  zm=  < 

^  Z2,  for  r  =  3,4  (  Z2,  for  9  =  3,4 

x,  for  r  =  1,3  (  x,  for  ?  =  1,3 

a  =  |  ,  (3=1 

y,  for  r  =  2,4  I  y,  for?  =  2,4 


(25a) 


(256) 


The  determinantal  equation  for  the  propagation  constant  kv o  can  be  solved  by  setting 
the  determinant  of  the  coefficient  matrix  of  (22)  equal  to  zero  : 


det[z(u/,fcvo)]  =0 
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B^Two  Symmetrical  Signal  Striplines  Loaded  with  Crossing  Strips 

In  this  section,  we  consider  the  case  where  two  identical  signal  striplines  of  width  u>i 
are  located  symmetrically  at  (±1^^)  as  shown  in  Fig. 3. 

For  the  even  modes,  Jmz(x,  y )  is  an  odd  function  of  x,  and  Jmy(x,y)  is  an  even  function 


of  x.  Therefore,  the  surface  currents  can  be  expanded  as 

N,  N, 

Mz,y)  =  x  Y  aJf[f{x,y)  +  y  ^  tyAffav)  (27a) 

j=-<V,  j=-Ni 

Ns  N« 

M*iy)  =  *  +  y'52dif(i'\x<y)  (276) 

;= l  ;=o 

where  the  basis  functions  are 

/J;)(*,y)=  [p,(*-Xc,u»i)  +  P1(*  +  *c,  tm)]eifc»*  (28a) 

/^e)(*,y)=  [To(*-a:c,t1;1)  +  ^o(*  +  *c,^l;1)]ei*«,'  (286) 

=  P;(x,Ic)(p/27r)T0(y,u;c)  (28c) 

/ii>(as»  y)  =  Tj{x,Le)Q(y,we).  (28  d) 


Following  the  same  procedure  as  for  the  case  of  one  signal  stripline  by  applying  the 
Galerkin’s  method,  a  determinants!  equation  similar  to  (26)  is  obtained. 

For  the  odd  modes,  Jmz(x,y)  is  an  even  function  of  x,  and  Jmy{z,y)  is  an  odd  function 
of  x.  Therefore,  the  surface  currents  are  expanded  as 

N ,  N, 

•7i(*>y)  =  *  Y  <*>/!?(*>»)  +  y  Y  bifi?(x'y)  (29a) 

j=-Nx  ;=-N, 

Nj  N4 

M*,y)  =  /!?(*,  y) +  y^rf;/i°)(*,y)  (296) 

j=l  j= 0 

where  the  basis  functions  are 


Pi(*  -  *c,wi)  -  Pi(*  +  *c,tnj)  e***J * 


(30a) 
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A;  (*1!/)—  —  *cj^l)  ~  Tq(z  +  ®e,t«l)j  e'**'* 

/&W)  =  ^(a?,Zc)(p/2ff)7’o(y^c) 

/l;)(*>y)  =  ^(x,Zc)Q(y,Wc) 


(306) 

(30c) 

(30<f) 


where 


m,v)  =  { 


Vj(t>v)  = 


f  (1/77)  cos[(2j  -  1  M/77],  -77/2  <  {  <  77/ 2 


0, 


elsewhere 


sin[(2;  -  l)irZ/rf]/y/(rf/2)»  -  £2,  -77/2  <  £  <  77/2 


0,  elsewhere. 

Let  Uj(kx,Tj)  and  Vj(kx,-q)  be  the  Fourier  transform  of  Uj(x,Tj)  and  V;(x,T7),  then 

U  (k  77 )  =  (2J  ~  1)(~1)J  cos(fet77/2)  _  - 
’  (*,>))’  -  [(2j  -  l)ir|*  _tf)< 


(31a) 

(316) 


(32a) 


^(*••1)  =  J- [-MM/2  -  (j  -  l/2)»)  -  J0(W2  +  0  -  1/2)*)]  = 


(326) 


When  <5*77/2  approaches  ±(j  —  l/2)7r,  Uj(k9,r 7)  approaches  l/(47r). 
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f!lQru»  Signal  Stripline  Loaded  with  Crossing  Strips  of  Infinite  Length 

In  this  section,  we  consider  the  case  where  one  signal  line  is  loaded  by  infinitely  long 
crossing  strips  as  shown  in  Fig. 2  with  Le  — »  oo.  When  the  crossing  strips  are  very  long 
such  that  the  reflections  from  ends  can  be  neglected,  we  can  assume  travelling  waves  along 
it.  So,  we  investigate  the  possibility  of  the  existence  of  such  a  mode  of  operation,  and  its 
effect  on  the  propagation  characteristics  of  the  signal  line. 


The  surface  current  on  the  signal  strip  is  of  the  same  form  as  in  the  case  (A)  of  finite 
crossing  strips.  For  the  crossing  strips,  we  choose  travelling  wave  basis  functions[9],  [10] 
and  some  local  basis  functions  on  the  center  to  account  for  the  effect  of  the  presence  of 
the  signal  line.  Hence,  the  surface  currents  are  expanded  as 

■V,  Nj 

Ji(z,y)  =  z  aifi?(x’y)  +  y  ^2  hfSjfav)  (33a) 

}=-N%  i=-N 1 

Mx,y)  =  +  y'%2djf(4?{z,y)  (336) 

i= i  ;=o 


where  the  basis  functions  f[^(z,y)  and  /j^(*,y)  are  the  same  as  fij{z,y)  and  f3j{x,y) 
respectively.  The  functional  forms  of  fl^{x,y)  and  f!^(z,y)  are 

I  [*,■(*,*)  -  «i(-*,*)](p/2ir)r,(»,»e),  1  <  j  <  Ni  -  1 

[-5m(Jfcez  -  w/2)  +  iSm(ktz)  (34a) 

+ Sm(-ktx  -  it  12)  -  i5m(-fce*)](p/2w)r0(y,u;c),  j  =  N3 


’  Ro{x,h)Q(y,we ),  j  =  0 

[ilj(*,fc)  +  Rj(-z,h)^Q(y, tne),  1  <  j  <  W4  -  1 

^-5m(fce*  -  w/2)  +  iSm{ktx) 

-Sm(-ktz  -  it/ 2)  +  iSm(-ktz)^Q(y,we)y  j  =  N4 


(346) 


where  ke  is  assumed  to  be  the  propagation  constant  of  a  single  crossing  strip  of  infinite 
length  in  the  absence  of  the  signal  line,  Rj(z,h)  is  the  local  basis  function  with  width 
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2h,  Sm(£)  is  the  travelling  wave  basis  function  with  m/2  periods  as  shown  in  Fig.4.  The 
explicit  forms  are 


sin  kt(h  —  \x  -  jh\)/  ainkeh,  ( j  -  1  )h  <  z  <  (j  +  1  )h 


0, 


elsewhere 


(35a) 


{sin£,  0  <  (  <  mir 

(356) 

0,  elsewhere. 

It  will  be  shown  that  only  finite  number  of  periods  of  the  travelling  wave  basis  functions 
are  sufficient  for  the  convergence  of  the  solution.  Any  increase  in  the  number  of  periods  of 
these  basis  functions  will  have  a  negligible  effect  on  the  numerical  results. 


Let  Rj(ka,w)  and  Sm(fci)  be  the  Fourier  transform  of  Rj(z,h )  and  Sm(£),  then 


R,(kt,h) 

Sm(kt) 


kee  lk*ik(cos  kzh  —  cos  keh) 


it  ain  kth(kl  —  k*) 


k* 


2 *(K  -  *2) 


[(-i) 


m  —i(m-K/kt)hm 
c 


-4 


(зба) 

(збб) 


When  kx  approaches  ±ke,  Rj(kZyh)  and  Sm(kx )  approach  (h/2n)e^lk*^K  and  2fam/4£e, 


respectively. 


With  these  basis  functions,  the  Fourier  transform  of  the  surface  current  Jm(kx,kyn)  can 
be  derived  as 

N,  N, 

•?.(*„*,.)  =i  £  6,,)  +  $  £  ‘if£’(*..*,»)  (37") 

j=-N  i  j——Ni 

Ns  Na 

Mk. ,  *,„)  =  i  £  CjFi?(k„k,„)  +  y  £  (376) 


J  =  1 


1=0 


where  F^(h*,fctn)  is  the  Fourier  transform  of  fjj\z,y),  i  =  1,2, 3, 4.  The  explicit  forms 
sure 


Fif(k„kyn)  =  SjnP1(katw ,)  (38a) 

Fif(kt1kyn)  =  SjMk^Wi)  (386) 

Fi?(kXik,n)  =  f„(*,n,u;e)A,>(M  (38c) 

Fi?(kmyk,n)  =  <j(*,n,toe)i4>(*,)  (38d) 
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where 


A.sj(ka) 


A*  j(kt) 


(  2R(;\k„h), 

1  kt[n(kl  -  A:*)]-1  [(-l)m* sin[fc«7r(m  +  1/2 )jk€) 
+(  —  l)m  sin(mirkx/ke)  —  i  sin(&,7r/2fce)j , 

{  R(o‘\k',h), 


1  <  j  <  Ni  -  1 


(39a) 

j  =  N> 
j  =  0 


2R('\kt,h),  1  <j<Nt 

kt[ir(kl  -  kl)]~l  j^-(-l)m  cos[A:I7r(m  4-  l/2)/fce] 

+t(-l)m  cos(mnkx/ke)  +  coa(kxir/2ke)  —  ij ,  j  =  jV4 


1 

(396) 


where  R^\kxth)  and  R^°\kx,t i)  are  the  even  and  the  odd  parts  of  Rj(kx,h )  respectively. 
It  is  observed  that  A3j(—kx)  =  —  .^(A:*),  and  ^4j(— kx)  =  A4j(kx). 


Applying  the  Galerkin’s  procedure  with  the  following  testing  functions  : 


=  Pi(*,w1)e~tfc«* 

(40  a) 

W7  j{x,y)  =  ^o(*,u>l)e~’*»>,, 

(406) 

™3j(*,y)  =  Rj(x,h)(p/2n)T0(y,wc) 

(40c) 

=  Rj(xi  h)Q(y,  we) 

(40<£) 

the  matrix  eigenvalue  equation  is  obtained. 
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IV-Numerical  Results  and  Discussions 

By  utilising  the  symmetry  properties  of  the  dyadic  Green’s  function,  the  basis  functions, 
and  the  testing  functions,  each  matrix  element  in  (22)  can  be  reduced  to  an  integral  over 
0  <  ka  <  oo.  In  computing  the  integrals  (23)  numerically,  the  path  of  integration  in  the 
complex  ka -plane  is  deformed  below  the  real  axis  to  avoid  the  poles  corresponding  to  the 
waveguide  modes[ll]. 

Fig.5  shows  the  dispersion  relation  for  a  single  signal  line  with  crossing  strips.  Numerical 
computations  were  performed  with  two  different  numbers  of  basis  functions,  and  the  results 
were  found  to  be  the  same  up  to  three  decimal  points.  The  basis  functions  used  are 
given  by  (15).  For  Le  =  2.3  mm,  the  first  stopband  occurs  in  the  frequency  range  when 
0.3162  <  kop/ir  <  0.3203. 

Fig. 6  shows  the  interaction  of  n  =  —1  Floquet  mode  with  TEi  (TMi)  parallel- plate 
waveguide  mode.  For  frequencies  above  /„,  ky  starts  to  have  large  imaginary  part,  giving 
rise  to  a  higher  order  stopband.  However,  we  are  interested  in  operating  frequencies  where 
ky  is  real  within  the  passbands  below  /,,  and  thus  the  region  above  /,  is  of  no  practical 
importance. 

Next,  the  effects  of  crossing  strip  length  Lc  on  the  lower  and  upper  ferquency  bounds 
of  the  stopband  are  investigated.  The  normalized  frequency  for  the  two  bounds  of  the 
stopband  is  presented  as  a  function  of  Le.  The  result  for  p  =  0.5  mm  is  plotted  in  Fig.7(a). 
It  is  observed  that  both  bounds  of  the  stopband  are  very  sensitive  to  the  crossing  strip 
length  Le •  This  behavior  is  repeated  when  Lc  is  changed  by  approximately  an  integral 
number  of  wavelengths.  This  can  be  explained  in  the  following  way  :  The  crossing  strips 
behave  like  open-circuited  stubs  periodically  loading  the  signal  line.  The  crossing  strips 
will  have  a  capacitive  or  inductive  behavior  depending  on  its  length.  At  a  certain  length  of 
crossing  strips,  the  behavior  switches  from  being  inductive  (  or  capacitive  )  to  capacitive 
(  or  inductive  ).  This  switching  occurs  at  Le  «  nAj,  where  is  the  wavelength  in  the 
dielectric  medium  calculated  at  the  center  frequency.  At  these  lengths,  the  stopbands 
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become  very  wide. 


In  Fig.7(b),  the  normalized  frequency  of  the  bounds  of  the  stopband  are  plotted  as 
function  of  Le  with  the  period  p  =  1.0  mm.  Behavior  similar  to  that  in  Fig. 7(a)  is 
observed,  but  the  values  of  Lc  at  which  the  switching  of  frequency  occurs  are  doubled. 

In  Fig. 8,  the  effects  of  the  crossing  strip  width  wc  on  the  stopband  frequency  bounds 
are  investigated.  The  normalized  frequency  at  the  bounds  of  the  stopband  is  presented 
for  Le  —  2.7  mm.  As  the  crossing  strip  width  becomes  smaller,  the  stopband  becomes 
narrower. 

In  Fig. 9,  we  investigate  the  effect  of  the  separations  <i  and  /j  on  the  stopband  while 
keeping  constant.  It  is  observed  that  the  stopband  becomes  smaller  when  the  separation 
is  decreased,  and  when  the  separation  is  larger  than  0.2  mm,  the  upper  frequency  bound 
of  the  stopband  reaches  a  constant. 

In  Fig.  10,  the  bounds  of  the  first  stopband  is  plotted  as  a  function  of  the  distance  t2 
while  fixing  the  separation  =  t3  =  constant.  It  is  observed  that  for  Le  —  1.0  mm,  the 
separation  t?  affects  the  upper  bound  of  the  first  stopband  significantly. 

Fig.  11  shows  the  case  of  two  coupled  signal  striplines  in  the  presence  of  periodic  crossing 
strips  of  finite  length.  The  frequency  bounds  of  stopband  are  presented  in  Fig.ll(a)  and 
Fig.ll(b)  for  the  even  mode  and  the  odd  mode,  respectively,  with  p  =  0.5  mm  and  Lc  =  1.7 
mm.  The  basis  functions  used  are  given  by  (27)  and  (29)  for  the  even  and  the  odd  modes, 
respectively.  When  the  separation  becomes  larger  than  Lc,  the  stopband  width  of  the  even 
mode  approaches  zero,  but  the  stopband  width  of  the  odd  mode  is  still  finite.  This  is 
because  the  odd  mode  has  stronger  coupling  between  two  signal  lines  than  the  even  mode. 

Fig.  12  shows  the  dispersion  relation  of  a  single  stripline  in  the  presence  of  crossing 
strips  of  infinite  length.  The  basis  functions  used  are  given  by  (33).  We  choose  the 
travelling  wave  basis  functions  to  have  three  periods.  The  results  using  seven  periods  are 
also  shown  for  comparison,  and  it  is  found  that  the  travelling  wave  basis  function  of  three 
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period*  is  sufficient.  The  imaginary  part  of  the  propagation  constant  is  approximately  a 
linear  function  of  frequency,  and  the  magnitude  can  be  as  high  as  one  percent  of  the  real 
part.  This  is  due  to  the  assumption  that  the  surface  current  along  the  crossing  strips  is 
a  travelling  wave.  Part  of  the  power  along  the  signal  line  couples  to  the  crossing  strips, 
exciting  travelling  wave  surface  current  flowing  away  from  the  signal  stripline.  These 
travelling  wave  surface  current  guides  some  power  away  from  the  signal  line,  hence  reduce 
the  guided  power  alour  the  signal  line. 

In  this  case,  the  passband  stopband  behavior  which  is  characteristic  of  periodic  struc¬ 
tures  does  not  appear.  The  wave  number  ky o  has  nonzero  imaginary  part  over  all  fre¬ 
quencies.  This  is  due  to  the  power  guided  by  the  travelling  wave  along  the  crossirg  strips. 
Around  ky 0  =  nir/p ,  the  separation  between  two  neighboring  crossing  strips  :s  nA/2  where 
A  is  the  wavelength  of  the  guided  mode.  The  power  carried  by  the  crossing  strips  at  these 
frequencies  is  very  small  because  the  current  on  the  signal  line  have  opposite  phase  on 
both  sides  of  the  crossing  strip.  Also,  in  this  case,  the  higher  order  waveguide  mode  is  not 
excited. 

Conclusions 

A  rigorous  dyadic  Green’s  function  formulation  for  the  periodic  structure  is  derived  to 
study  the  dispersion  properties  of  single  and  coupled  signal  lines  periodically  loaded  with 
crossing  strips.  The  passband  and  stopband  characteristics  are  investigated  when  crossing 
strips  are  of  finite  or  infinite  length. 

For  crossing  strips  of  finite  length,  the  stopband  properties  are  mainly  affected  by  the 
period,  the  length  of  crossing  strips,  and  the  separation  between  the  signal  and  crossing 
strips.  Also,  at  higher  frequencies,  higher  order  stopbands  occur.  For  crossing  strips  of 
infinite  length,  attenuation  along  the  signal  line  exists  over  the  whole  frequency  range  due 
to  the  power  guided  by  the  travelling  wave  along  crossing  strips. 


18 


References 


[1]  T.  Kitasawa  and  R.  Mittra,“An  investigation  of  striplines  and  fin  lines  with  periodic  stubs,” 
IEEE  Tran s.  Microwave  Theory  Tech.,  vol.MTT-32,  no. 7,  pp. 684-688,  July  1984. 

[2]  Y.  Fukuoka  and  T.  Itoh, “Slow- wave  coplanar  waveguide  on  periodically  doped  semiconduc¬ 
tor  substrate,”  IEEE  Trans.  Microwave  Theory  Tech.,  vol.MTT-31,  no.12,  pp. 1013-1017, 
December  1983. 

[3]  Q.  Gu  and  J.  A.  Kong,  “Transient  analysis  of  single  and  coupled  lines  with  capacitively 
loaded  junctions,”  IEEE  Trans.  Microwave  Theory  Tech.,  vol.MTT-34,  no. 9,  pp. 952-964, 
September  1986. 

[4]  B.  J.  Rubin, “The  propagation  characteristics  of  signal  lines  in  a  mesh-plane  environment,” 
IEEE  Trans.  Microwave  Theory  Tech.,  vol.MTT-32,  no. 5,  pp. 522-531,  May  1984. 

[5]  C.  H.  Chan  and  R.  Mittra,  “The  propagation  characteristics  of  signal  lines  embedded  in 
a  multilayered  structure  in  the  presence  of  a  periodically  perforated  ground  plane,”  IEEE 
Trans.  Microwave  Theory  Tech.,  vol.MTT-36,  no. 6,  pp. 968-975,  June  1988. 

[6]  N.  V.  Nair  and  A.  K.  Mallick,  “An  analysis  of  a  width-modulated  microstrip  periodic 
structure,”  IEEE  Trans.  Microwave  Theory  Tech.,  vol.MTT-32,  no.2,  pp. 200-204,  Febru¬ 
ary  1984. 

[7]  F.  J.  Glandorf  and  I.  Wolff,  “A  spectral-domain  analysis  of  periodically  nonuniform  mi¬ 
crostrip  lines,”  IEEE  Trans.  Microwave  Theory  Tech.,  vol.MTT-35,  no.3,  pp. 336-343, 
March  1987. 

[8]  J.  A.  Kong, u Electromagnetic  Wave  Theory ,”  New  York  :  Wiley,  1986. 

[9]  R.  W.  Jackson  and  D.  M.  Pozar,  ”Full-wave  analysis  of  microstrip  open-end  and  gap 
discontinuities,”  IEEE  Trans.  Microwave  Theory  Tech.,  vol.MTT-33,  no. 10,  pp. 1036-1042, 
October  1985. 

[10]  H.  Y.  Yang  and  N.  G.  Alexopoulos,  “A  dynamic  model  for  microstrip-slotline  transition  and 
related  structures,”  IEEE  Trans.  Microwave  Theory  Tech.,  vol.MTT-36,  no.2,  pp. 286-293, 
February  1988. 

[11]  E.  H.  Newman  and  D.  Forrai,  “Scattering  from  a  microstrip  patch,”  IEEE  Trans.  Antennas 
Propagat.,  vol.AP-35,  no.3,  pp. 245-251,  March  1987. 


Fig.l  Geometrical  configuration  of  signal  strip  lines 
periodically  loaded  with  crossing  metallic  strips 
embedded  in  layer  (/)  of  a  stratified  medium. 

Fig. 2  Geometrical  configuration  of  one  signal  strip  line 
periodically  loaded  with  crossing  strips 
embedded  in  a  one-layer  medium. 

Fig.3  Geometrical  configuration  of  two  signal  strip  lines 
periodically  loaded  with  crossing  strips 
embedded  in  a  one-layer  medium. 

Fig.4  The  basis  functions  used  for  infinitely  long  crossing  strips, 
m  =  4,  iVj  =  N*  =  6. 

Fig. 5  Dispersion  relation  of  one  signal  strip  line  periodically 
loaded  with  crossing  strips  of  finite  length, 
er  =  10,  <1  =  t2  =  <3  =  0.2  mm,  p  =  0.5  mm, 
wi  =  we  —  0.125  mm,  Lc  =  2.3  mm, 

Nx  =  1,  N3  =  0,  N3  =  4,  N<  =  3, 
x  :  Ni  =  1,  N2  =  0,  N3  =  6,  N*  =  5, 
u>u  :  upper  bound  of  the  first  stopband, 
uji  :  lower  bound  of  the  first  stopband. 

Fig.6  Interaction  of  Floquet  modes  with  TEi  (TMj) 
parallel-plate  waveguide  mode, 

k  =  <jJ yJlXQlQtr,  d  =  <1  +  tj  +  tj. 


Fig.7(a)  The  effects  of  Le  on  the  upper  and 
lower  bounds  of  the  stopband, 

€r  =  10,  <1  =  t2  =  t3  =  0.2  mm, 
p  =  0.5  mm,  ujj  =  we  —  0.125  mm, 
Ni  =  1,  N?  =  0,  N3  -  5,  N<  =  4.  . 

Fig. 7(b)  The  effects  of  Le  on  the  upper  and 
lower  bounds  of  the  stopband, 
er  =  10,  ti  =  *2  =  t3  =  0.2  mm, 
p  =  1.0  mm,  w3  =  we  =  0.125  mm, 
N3  =  1,  N2  =  0,  iV3  =  5,  N4  =  4. 
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Fig. 8  The  effects  of  toc  on  the  upper  and 
lower  bounds  of  the  stopband, 

«r  =  10,  <1  =  =  *s  =  0.2  mm, 

p  =  0.5  mm,  tui  =  0.125  mm,  Le  =  2.7  mm, 
Nj  =  1,  N2  =  0,  .V3  =  4,  jV4  =  3. 


Fig.9  The  effects  of  and  t3  on  the  upper  and 
lower  bounds  of  the  stopband, 
er  —  10,  =  0.2  mm,  p  =  0.5  mm, 

w\  =  toe  =  0.125  mm,  Lc  =  1.0  mm, 

JV3  =  1,  N,  =  0,  N3  =2 ,N<  =  1. 


Fig.10  The  effects  of  t3  on  the  upper  and 
lower  bounds  of  the  stopband, 

€r  =  10,  =  f3  =  0.2  mm,  p  =  0.5  mm, 

w\  =  wc  ss  0.125  mm,  £e  =  1.0  mm, 

JVi  =  1,  JV3  =  o,  N3  =  2,  iV4  =  1. 


Fig.ll(a)  The  effects  of  xc  on  the  upper  and 
lower  bounds  of  the  stopband, 
for  two  signal  strip  lines, 
er  =  10,  t\  =  t3  —  t3  —  0.2  mm,  p  =  0.5  mm, 
<03  =  106  =  0.125  mm,  Lc  =  1.7  mm, 

N\  =  1,  N2  =  0,  N3  =  3,  JV4  =  2. 


Fig.ll(b)  The  effects  of  xe  on  the  upper  and 
lower  bounds  of  the  stopband, 
for  two  signal  strip  lines, 
er  =  10,  <i  =  <3  =  t3  =  0.2  mm,  p  =  0.5  mm, 
w\  =  Wc  =  0.125  mm,  Lc  —  1.7  mm, 

N3  =  1,  Nt  =  0,  iVj  =  3,  Na  =  2. 


Fig. 12  The  dispersion  relation  for  one  signal  strip  line 

periodically  loaded  with  crossing  strips  of  infinite  length, 
er  =  10,  tj  =  f3  =  t j  =  0.2  mm,  p  =  2.5  mm, 
u>j  =  toe  =  0.125  mm,  h  =  A,/8,  m  =  6, 
x  :  m  =  14. 
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Abstract-  A  new  method  for  analyzing  frequency-dependent  transmission  line  systems 
with  nonlinear  terminations  is  presented.  The  generalized  scattering  matrix  formulation  is 
used  as  the  foundation  for  the  time  domain  iteration  scheme.  Compared  to  the  admittance 
matrix  approach  proposed  in  a  previous  paper,  it  has  the  advantage  of  shorter  impulse 
response  which  leads  to  smaller  computer  memory  requirement  and  faster  computation 
time.  Examples  of  a  microstrip  line  loaded  with  nonlinear  elements  are  given  to  illustrate 
the  efficiency  of  this  method. 

I.  INTRODUCTION 

In  recent  years,  the  effect  of  the  interconnection  lines  on  high-speed  integrated 
circuits  has  become  more  and  more  important.  As  the  speeds  of  integrated  circuits 
increase,  the  propagation  delay  as  well  as  the  dispersion  and  loss  of  interconnection 
lines  can  no  longer  be  neglected.  Traditional  lumped  element  circuit  models  must 
be  supplemented  by  dispersive  transmission  line  models  in  order  to  account  for 
these  effects.  This  has  created  the  need  for  new  numerical  procedures  that  can  be 
easily  incorporated  into  current  CAD  tools.  To  make  matters  more  complicated, 
the  interconnection  lines  are  terminated  with  not  only  linear  elements  but  also 
nonlinear  semiconductor  devices,  such  as  diodes  and  transistors. 

Several  techniques  are  now  commonly  used  to  deal  with  nonlinear  circuit  prob¬ 
lems,  for  example,  the  direct  time  domain  approaches  [1,2],  and  the  semi-frequency 
domain  approaches,  such  as  the  harmonic  balance  [3,4]  and  the  modified  har¬ 
monic  balance  techniques  [5,6j.  Semi-frequency  domain  approaches  are  useful 
for  microwave  and  millimeter  wave  integrated  circuits  but  become  impractical 
for  digital  integrated  circuits  because  of  their  wide-band  nature.  On  the  other 
hand,  frequency-dependent  parameters  often  make  it  awkward  to  apply  the  direct 
time  domain  approach  to  interconnection  line  systems.  The  time-domain  finite- 
difference  method  [7]  and  the  time-domain  method  of  moments  [8]  have  been 
proposed  to  deal  directly  with  electromagnetic  scattering  from  nonlinear  loads. 
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However,  dispersion  problems  are  absent  from  the  discussions. 

Liu  and  Tesche  [1]  developed  a  combined  time-domain  frequency-domain  treat¬ 
ment  of  antennas  and  scatterers  with  nonlinear  loads.  In  their  work,  the  transfer 
function  (impulse  response)  of  the  linear  portion  of  the  investigated  system  is  first 
obtained  through  the  frequency- domain  analysis,  and  it  is  then  used  to  solve  the 
entire  nonlinear  problem  in  the  time  domain.  Subsequent  improvements  to  this 
method  have  been  suggested  by  Djordjevic,  Sarkar  and  Harrington  [9]  through 
artificially  introducing  pairs  of  quasi-matched  passive  networks,  and  by  Canig- 
gia  [10]  and  Schutt-Aine  and  Mittra  [11],  through  macromodel  and  scattering 
parameter  analysis  based  on  a  fixed  reference  impedance. 

In  this  paper  we  shall  present  an  extended  and  more  natural  method  that 
will  completely  eliminate  the  need  for  any  artificial  networks  or  fixed  reference 
impedances.  Its  close  ties  to  the  physics  of  transmitted  and  reflected  waves  on 
transmission  lines  also  help  in  achieving  the  purpose  of  reducing  computation 
time.  The  algorithm  is  explained  in  the  next  section.  In  Section  III  the  de¬ 
tails  of  applying  our  formulations  to  frequency-dependent  transmission  lines  with 
nonlinear  loads  are  illustrated  in  the  analyses  of  a  nonlinearly-loaded  dispersive 
transmission  line. 

n.  ANALYSIS  BASED  ON  WAVE  TRANSMISSION  AND  REFLECTION 

An  arbitrary  system  of  n  dispersive  transmission  lines  can  be  represented  by  the 
following  coupled  linear  ordinary  differential  equations  in  the  frequency 

-slv'l=M£i'w  +  1*,'[/l 

-^W=MC]-[v]  +  |g]-(vi 

Treating  the  n  transmission  lines  as  a  2n  -port  system,  we  can  derive 
the  admittance  matrix  [K] ,  which  relates  terminal  voltages  to  terminal 

2n 

h  =  £  W  <2> 

4=1 

The  time  domain  counterparts  become  convolution  relations: 

*;(*)  =  ]l/  dryjk(i-r)v 4(r)  (3) 

4=1 Jo 

where  yjj,  is  the  inverse  Fourier  transform  of  The  terminal  voltages 

and  currents  for  any  particular  system  can  then  be  uniquely  determined  once  the 
terminal  conditions 

M‘)l  =  [/("(<))]  (4) 

are  given.  If  all  the  terminations  contain  only  linear  elements,  we  can  solve  the 
problem  in  the  frequency  domain.  Otherwise,  iteration  procedures  are  usually 
required  for  obtaining  the  solutions.  The  analyses  presented  in  [1]  and  [9]  are  based 
on  equations  (3)  and  (4).  Although  their  approaches  are  straightforward,  there 


domain: 

(1) 

from  (1) 
currents: 


2Vaw«tenf  iM/|<i<  of  Tmnsmiuion  Line  Salterns 
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exist  problems  that  could  possibly  affect  the  efficiency  of  numerical  computation. 
First  of  all,  the  parameter  yyfc(f)  is  equivalent  to  the  current  measured  at  port  j 
with  a  voltage  impulse  excitation  of  unit  amplitude  at  port  k  while  ports  other 
than  k  are  short-circuited.  Owing  to  strong  reflections  at  the  terminations,  the 
duration  of  yy*( <)  is  usually  long  for  slightly  lossy  transmission  line  systems, 
and  even  infinite  for  a  lossless  system.  The  long  duration  puts  great  demands  on 
computer  memory  and  execution  time  as  required  to  perform  convolution  integrals 
of  [v(f)]  and  [y(f)j  during  the  iterative  solution  of  nonlinear  equations. 

In  order  to  overcome  the  disadvantage  of  using  the  parameters  {yy*},  the  au¬ 
thors  in  [9]  artificially  insert  a  pair  of  complementary  passive  networks  between 
the  end  of  the  transmission  lines  and  the  actual  terminations.  The  purpose  of 
that  is  to  make  the  augmented  linear  network,  which  consists  of  the  transmission 
line  portion  and  the  artificial  network  directly  connected  to  the  end  of  the  trans¬ 
mission  lines,  a  quasi-matched  linear  system  so  that  the  duration  of  the  impulse 
responses  for  the  augmented  linear  system  can  be  effectively  shortened.  However, 
the  other  artificial  network  which  is  directly  attached  to  the  original  load  will 
contain  negative  resistors  and  hence  may  render  the  numerical  solution  unstable, 
especially  when  the  transmission  lines  are  lossless. 

Instead  of  dealing  with  terminal  voltage  and  current,  we  will  analyze  the  trans¬ 
mission  line  system  from  the  viewpoint  of  voltage  waves.  We  choose  the  input 
and  output  waves  at  the  terminal  ports  of  the  transmission  lines  as  the  variables 
of  the  problem  as  shown  in  Fig.  1.  The  parameters  {By}  and  {Cy}  are  defined 
as  follows: 


(5) 

c,H  [v,W  +  z0iH .  /,(*)] 

(6) 

Frequency-Dependent 


Transmission  Line  System 


c2 

b2 


C2i 


C2n 

b2n 


Figure  1.  Linear  multi-port  network  consisting  of  transmission  lines. 
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where  Zoj(u)  =  yJLjj{.^)/Cjj(uj)  is  the  frequency-dependent  characteristic  impe¬ 
dance  on  line  connected  to  port  j  .  The  linear  dispersive  transmission  line  system 
is  thus  characterised  by  a  scattering  matrix  [5jj] ,  i.e., 


r  Ci  i 

5n 

S 12 

Sn 

*  •  *  5i>2n  ' 

c2 

■$21 

S22 

^23 

•  •  •  52l2n 

Ci 

= 

^31 

s22 

533 

•  •  •  Sit2n 

.c2n. 

.52n,i 

S2n,2 

52n,3 

**•  52n,2». 

B  i 
B2 

lB2nl 


(7) 


or 


2n 


(*  =  1,2,..  .,2n) 


(8) 


Cj{w)  =  ^2  Sjk(ui)Bk(u;) 

k= 1 

It  is  easy  to  realize  from  the  above  equations  that  for  all  j  k ,  Sjfc(u>)  is 
equal  to  2Vj(u>)  if  all  ports  axe  loaded  with  their  transmission  line  characteristic 
impedances  and  only  a  voltage  source  of  amplitude  1  is  applied  at  port  k .  This 
corresponds  to  impulse  response  or  transfer  functions  in  the  time  domain.  If  the 
coupling  between  individual  transmission  lines  is  weak,  the  system  will  be  close 
to  being  perfectly  matched.  In  this  case  we  can  conclude  that  the  inverse  Fourier 
transform  of  Sj k  ,  denoted  as  hjk(t) ,  will  be  of  much  shorter  duration  than  yn(t) . 
Therefore  we  can  effectively  reduce  the  memory  required  to  store  the  past  values 
of  hjt  and  the  time  to  compute  the  convolution  integrals  in 

cj(<)  =  J2  /  dThjk(i-T)bk(r)  (9) 

*=1  J0 

without  inserting  any  artificial  networks. 

We  now  have  to  solve  [b(t)\  and  (c(l)j  by  incorporating  the  nonlinear  boundary 
conditions  of  the  problem.  Specifically  at  port  j ,  equations  (5)  and  (6)  lead  to 


[8] 


and 


V>(«)  =  B,(»)  +  c,M 


(10) 


2C,(«)  =  ViH  +  Zo>(«)/J(w)  (11) 

By  taking  inverse  Fourier  transform  on  both  sides  of  (10)  and  (11),  we  obtain 
their  time  domain  counterparts: 

vji*)  =bj(t)  +  cj(t)  (12) 

Vj(i)  =  2 cj{t)  -  [  dr  z0j{t  -  r)  fj{vj{r))  (13) 

Jo 

where 

*0,(1)  =  [ZoiM] 

Our  problem  will  be  solved  in  a  time-marching  fashion.  At  any  instant  t ,  the 
iteration  procedure  is  as  follows: 

(i)  Set  up  initial  guess  of  (6(t)] .  A  reasonable  choice  is  to  take  values  from  the 
previous  time  step. 
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(ii)  Compute  individual  {cj(t)}  using  (9). 

(iii)  Apply  standard  nonlinear  equation  techniques  such  as  Newton-  Raphson  me¬ 
thod  to  (13)  to  solve  for  individual  {fj(t)}  ■ 

(iv)  Obtain  the  next  guess  of  [&(*)]  from  the  relation  bj(t)  =  vj(t)  —  cj(t),  and 
compare  the  guess  with  the  previous  one.  If  the  error  is  above  a  pre-set 
tolerance,  repeat  steps  (ii)  to  (iv)  with  the  new  guess. 

For  all  practical  purposes,  we  further  divide  these  variables  into  two  sets.  The 
first  corresponds  to  the  source  side  of  the  transmission  line  system  (Fig.  2),  iden¬ 
tified  by  odd-numbered  subscripts,  and  the  other  corresponds  to  the  load  side 
(sourceless)  with  even-numbered  subscripts  (Fig.  3).  The  reason  for  doing  so  is 
based  the  following  observations.  The  transfer  functions  linking  {62;}  to  {cjy+1} 
and  {bjj+i}  to  {c2j}  include  time  delay  operators  in  order  to  account  for  the 
finite  speed  of  propagation.  In  other  words,  {c2y+i}  depend  not  on  the  present 
but  the  past  values  of  {&2j}  and  the  like  for  {c2j}  on  {&2;+i}*  As  we  carry 
out  the  iteration  procedure  step  by  step,  the  present  values  in  one  set  will  not 
interfere  with  those  in  the  other.  Therefore,  the  update  of  variables  can  be  done 
simultaneously  for  any  given  time  if  parallel  processing  facilities  are  available. 

In  the  next  section  we  shall  present  the  application  of  this  method  to  a  single 
dispersive  line  loaded  by  nonlinear  impedance. 


m.  NONLINEARLY-LOADED  DISPERSIVE  TRANSMISSION  LINE 

Shown  in  Fig.  4  is  a  uniform  dispersive  transmission  line  of  length  l  driven  by 
a  source  eo(t)  with  a  linear  source  resistance  R$  at  one  end,  and  terminated 
by  a  nonlinear  resistor  An  at  the  other  end.  The  transmission  line  portion  of 
this  problem  can  be  described  in  terms  of  the  frequency  dependent  characteristic 
impedance  Zq(u>)  and  effective  propagation  constant /?(u/) .  It  can  be  shown  that 
the  frequency-domain  scattering  matrix  is  given  as  follows 


> 


188 


G%  et  al. 


"2  x 


Figure  3.  Even-numbered  ports  (load  side). 
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/2(0 
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Figure  4.  Nonlinearly  loaded  dispersive  transmission  line. 


’Sn(u/)  5i2(w)' 

0  e-i0(«)l  ‘ 

.52i(w)  522(u;), 

e-j0(v)l  o 

(14) 


Noting  that  both  5n  and  S22  are  zero  for  all  time;  therefore,  we  only  need 
to  consider  5i2  and  Sji  ,  which  involve  time  delays.  The  simplicity  of  (14) 
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can  be  attributed  to  our  definition  of  scattering  matrix.  As  an  example,  we 
consider  a  microstrip  line  of  width  W .  The  substrate  is  of  thickness  h  and 
has  a  dielectric  constant  tr .  Numerous  empirical  formulas  are  available  from 
the  literature  (see  for  example  (13-15J).  In  this  paper,  we  will  use  the  following 
expressions  to  calculate  the  frequency  dependence  of  microstrip  line  characteristic 
impedance  and  effective  dielectric  constant  ee(u>) ,  which  is  defined  from  /3(u/)  = 
Vy/t* 0<«<0  : 


€«(u/)  =  tr  ~ 


tT  —  e«(0) 


Zo{u)  = 


1  + 
r]h 


c«(0) 


(- 

\ut 


and  the  effective  width  We(u)  is  governed  by  the  equation 

Wn(Q)  -  W 


We(w)  =  W  + 


1  +  (Jj/ug 


where 


ttZo(O) 

Mo* 


(15) 

(16) 


(17) 

(18) 


“'9  Wsfc 

and  the  low-frequency  limits  of  ee  and  Zq  are  respectively 

«.(0)  =  ^  +  ~F{W/h) 

with 

P(W/h\  -  /  (1  +  12 h/W)-1'2  4-  0.04(1  -  W/h)2,  if  W/h  <  1 
K  7  ]  \  (1  +  12h/W)-V2,  if  W/h  >  1 

and 


(19) 

(20) 
(21) 


Zo(0)  =  { 


2*VuW) 


if  W/h  <  1 


1 7^  + 1,393 + 0-667U(t  + 1444)} 


-l 


if  W/h  >  1 


(22) 


Note  that  H^«(0)  in  (17)  is  derived  from  (16)  by  substituting  in  Zo(0)  and  e«(0). 

In  this  paper,  we  assume  that  the  width  W  is  equal  to  0.508  mm  and  the  depth 
and  the  dielectric  constant  of  the  substrate  are  equal  to  0.216  mm  and  10.2  respec¬ 
tively.  It  can  be  readily  derived  from  the  formulas  that  the  effective  permittivity 
in  the  low-frequency  limit  is  equal  to  7.46,  and  the  high-frequency  characteristic 
impedance  is  Z<j(oo)  =  50.2  Q .  The  calculation  of  the  impulse  response  function 
hl2  (t)  (=  hi i) ,  the  inverse  Fourier  transform  of  5u(u>) ,  however,  is  not  trivial. 

Because  lim  5ia(u>)  does  not  approach  0 ,  the  inverse  Fourier  transform  is  sin- 

v  >oo 

gular.  If  numerical  computations  are  not  carried  out  with  great  care,  the  accuracy 
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would  be  questionable.  In  [16]  a  similar  Fourier  integral  was  calculated.  The  only 
difference  is  that  the  input  is  not  an  impulse  function,  and  hence  it  is  numerically 
integrable.  The  authors  proposed  using  Taylor  expansion  method  and  the  method 
of  stationary  phase  for  narrow- band  input  signals.  But  for  other  input  waveforms 
a  brute  force  numerical  integration  method  was  applied  which  proved  to  be  very 
time-consuming.  To  overcome  the  difficulty  of  numerical  integration  once  and  for 
all,  our  approach  is  to  separate  the  transfer  function  $12  into  two  parts,  one  that 
can  be  analytically  inverted  to  an  impulse  function,  and  one  that  is  well-behaved 
and  integrable.  The  latter  requires  that  the  integrand  approaches  0  as  u  goes 
to  infinity.  This  leads  to  a  natural  way  of  separation: 

=  (e-jd(w)i  _  e~J0oo I'j  + 

-  e-)0aol  |  f  e-j(/3(w)-/3oo)f  _  l  j  +  1  }  (23) 

where  /?oo  =  <*>  • 

In  (23),  the  factor  corresponds  to  a  time  delay  of  rr  =  l/(cy/<c+)  and 

the  inverse  Fourier  transform  of  1  is  an  impulse  function.  Therefore,  defining 
h(t)  =  <^-1{e~j(/3(<‘,)-&»)i  —  1} ,  we  then  have 


hl2{t)  =  6{t  -  rr)  +  h(t  -  rr)  (24) 

Now  only  h(t)  needs  to  be  evaluated  numerically.  Because  the  kernel  of  this 
Fourier  inversion  has  essentially  a  finite  range  of  integration,  the  calculation  be¬ 
comes  easier.  The  impulse  response  function  for  a  10  mm  long  microstrip  line  is 
shown  in  Fig.  5.  Only  a  limited  portion  of  hi2(t)  surrounds  the  impulse  function. 
This  is  consistent  with  our  claim  that  the  impulse  response  has  a  very  limited 
duration. 

Because  the  scattering  parameters  h\i  and  h22  are  zero,  the  input-output 
wave  variable  pair  {6j(f),  Cj(t)}  are  only  linked  through  (12).  The  iterative  solu¬ 
tion  to  a  single  transmission  line  problem  is  therefore  relatively  simple.  Once  we 
finish  calculating  vj(t) ,  iteration  step  (iv)  will  readily  return  the  correct  values 
for  bj(t) .  There  is  no  need  to  go  back  to  step  (i). 

We  first  examine  the  response  from  a  Gaussian  pulse  with  an  amplitude  of  1.0 
volts  and  a  width  of  10  ps  measured  at  its  half  amplitude  is  used  a*  the  source. 
The  source  resistance  R$  is  50  fl .  The  load  characteristics  is  fully  described  by 
the  following  equation: 


i2  =  70  [«p  (0^5)  -  l]  nA 


(25) 
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Figure  5.  Impulse  response  J»i2(<)  • 
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Figure  6.  Terminal  voltages  vs.  time  at  both  the  source  end  and  the  load 
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Figure  0.  Load  voltages  for  different  line  lengths. 
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Figure  10.  Load  voltages  for  different  pulse  widths. 
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Figure  11.  Load  voltages  for  another  type  of  nonlinear  termination. 

For  the  case  Tq  =  1.0,  the  time  responses  at  both  ends  are  plotted  in  Fig.  6. 
The  pulse  has  been  broadened  and  the  negative  trailing  edge  is  rather  significant. 
In  order  to  analyze  what  comes  from  dispersion  of  the  transmission  line  and  what 
comes  from  the  load  nonlinearity,  we  compared  the  plots  of  the  load  responses 
when  Iq  is  varied  from  0.2  to  50.0  in  (25),  and  when  the  load  is  a  50  Cl  linear 
resistor.  The  results  are  depicted  in  Figs.  7  and  8.  Apparently  the  nonlinearity 
contributes  most  to  the  broadening  effect.-  The  negative  trailing  edge  originates 
from  dispersion,  but  is  enhanced  by  the  nonlinear  load.  Since  the  load  charac¬ 
teristics  is  similar  to  that  of  a  typical  diode,  it  behaves  like  an  open  circuit  with 
respect  to  an  incoming  negative  voltage  wave.  This  gives  rise  to  roughly  twice  the 
response  compared  to  the  one  at  the  50  Cl  load,  which  is  nearly  matched  to  the 
transmission  line.  The  large  negative  trailing  edge  also  accounts  for  the  extended 
ringing  after  it  is  reflected  back  to  the  source  end. 

In  Fig.  9,  the  load  end  voltages  for  5  mm  and  20  mm  lines  are  compared  against 
the  10  mm  line  case.  The  number  of  zero  crossings  increases  with  length,  as 
expected  from  the  transmission  line  dispersion.  Similar  phenomenon  is  observed 
when  shorter  pulses  are  injected,  as  in  Fig.  10.  The  centers  of  these  pulses  are 
intentionally  separated  to  allow  clearer  comparison.  It  is  interesting  to  note  the 
similarity  of  the  5  ps  pulse  output  and  the  20  mm  line  output  in  Fig.  9. 

We  have  also  tested  our  iteration  algorithm  on  another  class  of  nonlinear  loads 
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with  i-v  characteristics  described  as 


(26) 


Unlike  the  one  described  in  (25),  there  is  cap  on  the  positive  current.  The  result  is 
that  positive  voltage  has  a  larger  amplitude,  which  is  controlled  by  Iq  ,  as  shown 
in  Fig.  11. 

The  time  resolution  used  in  all  but  the  5  psec  pulse  input  case  is  0.4  psec  with 
a  total  of  1000  points.  For  the  latter,  the  resolution  is  0.2  psec  and  there  are  2000 
points.  Because  very  limited  portion  of  the  impulse  response  /112(f)  is  significant, 
the  actual  number  of  points  involved  in  the  convolution  integral  is  far  lower.  We 
used  the  Newton- Raphson  iteration  procedure  for  the  nonlinear  equation  (13). 
On  a  VAX  Station  3500,  the  testing  cases  take  about  4  to  27  seconds  of  CPU  time 
to  generate  the  solutions  with  1000  points.  The  large  variation  in  computation 
times  is  due  to  different  convergent  rates  for  different  types  of  nonlinearity. 


IV.  CONCLUSIONS 

A  new  method  for  the  transient  analysis  of  a  frequency-dependent  transmission 
line  system  terminated  with  nonlinear  loads  has  been  presented.  This  method  is 
not  only  effective  for  saving  the  CPU  time  required  for  solving  nonlinear  transient 
problems,  but  is  also  compact  and  natural  in  form.  Our  generalized  scattering 
matrix  approach  is  closely  tied  to  the  concept  of  waves.  Therefore,  no  extended 
reflection  arises  as  a  result  of  artificial  boundary  conditions  as  can  occur  with  the 
admittance  matrix  method,  and  duration  of  the  impulse  responses  for  the  waves 
in  the  transmission  line  system  is  very  limited.  This  is  the  key  to  reducing  the 
amount  of  computation  time  and  memory. 

The  detailed  procedure  for  solving  this  kind  of  nonlinear  transient  problem 
is  given  through  an  analysis  of  a  nonlinearly-loaded  microstrip  transmission  line 
with  linear  source  resistance.  Extension  of  this  approach  to  multiple  transmission 
line  systems  with  nonlinear  source  and  terminations  is  being  studied  and  will  be 
reported  in  future  work. 
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